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Introduction and Summary 

This report summarizes progress for the period of 
October 1, 1970 to March 31, 1971 performed under NASA 
Grant NGL 22-009-012, "Optical and Infrared Masers". Seven 
topics are summarized in the areas of high resolution studies 
of atoms and molecules, molecular excitation and relaxation 
mechanisms, short pulse experiments, and laser harmonic fre- 
quency mixing. Four appendices containing more detailed 
discussions are also included. 
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I . Application of Laser Induced Line Narrowing Techniques 

21 

to Measurement of Hyperfine Structure of Ne 

The laser-induced line narrowing technique has been used 

21 

to measure hyperfine structure in Ne . This technique, which 

utilizes laser-induced changes in the atomic state velocity 

distribution, permits the resolution of structure which is 

normally obscured by Doppler broadening. Measurements of the 

* 

narrowed spectral lines have yielded values for the isotope 

on oi on ° 

shifts of Ne , Ne , and Ne at 1.15u and 609 6A and hyper- 

21 

fine constants of Ne , which lead to an accurate value of the 
21 

Ne quadrupole moment. These measurements are described m 
detail in Appendix A. 
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XI . Excitation and Relaxation Mechanisms in the HF Laser 

Molecular Relaxation Studies in HF Gas 

Work has begun on assembling the necessary experi- 
mental apparatus for investigating vibrational and rotational 
relaxation processes occurring in pure, gaseous HF and in mix- 
tures of HF and other non-reactive gases. At present, a high 
vacuum, HF handling system has been designed and assembled. 

The system is constructed entirely of Kel-F and Teflon plastics. 
These materials are inert to chemical attack by HF, and, there- 
fore, the formation of impurities via chemical reaction is pre- 
vented. The system includes two cold traps for use in further 
purifying HF by successive distillation. Purity of the condensed 
sample can be ascertained by electrical conductivity measurements. 
Finally, a capacitance manometer, for precision pressure measure- 
ments was designed and fabricated. The manometer design is similar 
to that of a standard capacitance manometer, except that only non- 
reactive teflon materials were used for the interior of the sensing 
head. 


Near Infrared HF Laser 

A systematic investigation of the high energy, high gain 
pulsed HF laser is in progress. Initial experiments with a trans- 
verse discharge (TEA) HF, DF laser are described in a paper pub- 
lished in Applied Physics Letters and included as Appendix B of 
this report. Briefly, the laser is capable of high power, high 
energy (> 50 m.j.) pulsed output in the attractive 3 and 4u regions. 
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The laser device is simple and compact, and utilizes relatively 

t 

nontoxic gases. However, its spatially non-uniform excitation 
scheme and its high gain make single transverse mode operation 
difficult. To overcome these difficulties a short narrow bore, 
high pressure, longitudinal discharge laser has been designed 
and operated. It is capable of pulse energies of several milli 
joules and peak powers of >.2kW. Preliminary experiments indicate 
its single mode behavior is superior to that of the transverse dis- 
charge laser. A comparison of the spectral characteristic of the 
two lasers is now in progress. 

Long Wavelength HF Laser 

A wide bore, longitudinal discharge HF laser, to be used 
in studying far infrared emission in HF, has been built and operated. 
Laser excitation is accomplished by electrically pulsing a flowing 
mixture of and a fluorinated gas species such as SFg , Freon 14, or 
Freon 13. The device design is such that a wide range of laser opera- 
ting conditions and gas mixtures can be investigated. In studying 
the laser output spectrum, a fore-prism-grating spectrograph combina- 
tion has been found to be a convenient method of isolating near and 
far infrared laser lines. 

Using this apparatus, we have observed intense, pure rota- 
tional lasing in all three of the above named fluorinated gases. 

Laser wavelength ranged from 9 to 15 y; corresponding to rotational 
transitions from J=17 to J=30. Many of these laser wavelengths fall 
within the well-known 8 to 14 y atmospheric window. Studies of the 
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variation in laser spectral properties with gas pressure and 
gas mixtures have begun. Initial results have already provided 
clues to the process leading to rotational population inversion. 
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Ill . Coherent Effects in Short Pulse Propagation 


The study of fluorescence induced by short pulses of radia- 
tion in CC >2 has been concluded and the results submitted for pub- 
lication (see Appendix c). It has been demonstrated that one can 
observe the coherent excitation of an optical level in this way, 
even in a very weakly absorbing transition, and extract values of 
the matrix element and relaxation times. 

Present efforts in the area of short pulse propagation and 
coherent interaction have shifted to an investigation of adiabatic 
rapid passage. ^ This technique has long been used in magnetic 
resonance to invert a spin system and study relaxation processes. 
•Previous works in our laboratory on Stark spectroscopy of NH^ and 
NI^D has uncovered a number of transitions which can be swept through 
resonance with a C0 2 or ^0 laser transition. Because of their 
large absorption coefficient and tunabilitv, these transitions 
make excellent candidates for the study of adiabatic rapid passage. 
These experiments are of interest both from the point of view of 
studying the propagation of chirped laser pulses through an absorbing 
or amplifying medium, and for the study of relaxation processes. 

The laser and Stark cell are functioning properly and adequate 
absorption signals have been observed to make an attempt at obser- 
vation of adiabatic rapid passage feasible in the near future. 
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E. B. 'preacy, Phys. Letters 27A , 421 (1968). 



IV. Extension of Laser Harmonic Frequency Mixing Techniques . 

Efforts at frequency measurement in the infrared are 
currently directed at establishing a high accuracy laser fre- 
quency multiplying chain. The objectives are to, produce a pre- 
cise determination of the speed of light, superheterodyne 
detection at infrared wavelengths , and an accurate frequency 
measurement of the saturation resonances in CO 2 with possible 
application as a frequency standard. 

Based on the earlier work of Freed and Javan'*' , consider- 
able progress has been made (see VI below) in the observation 
of extremely narrow saturation resonances in C0 2 . The sharp 
resonances show considerable promise as high accuracy frequency 
standards. The laser frequency multiplier chain will be used 
for high accuracy frequency measurements of those resonances. 

A major obstacle to highly accurate frequency measurement 
with laser chains has been the difficulty of comparing several 
drifting laser frequencies to an absolute frequency standard. 

In 1968, Hocker, Small and Javan demonstrated in a point contact 
mixer that microwave frequency sidebands may be added to a laser 
carrier frequency and that drift in the laser may be offset by 
electronically tuning the microwave source. Work is now pro- 
ceeding toward construction of a laser frequency chain of elec- 
tronically stabilized sidebands which can be frequency locked 
to a frequency standard of absolute calibration. 

1 C. Freed & A. Javan, App. Phys. Letts. 17, 53 (1970). 

2 Quarterly Progress "Report, January ,22, 1969. 
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Figure I shows the method of producing a 9.3p stable fre- 
quency source which is phase-locked to a conventional frequency 
standard. A V-band microwave klystron, K^, is phase-locked to 
the frequency standard through a commercially' available frequency 
multiplier chain and phase-locked feedback loop, indicated as a 
lock box. The 12th harmonic of the klystron is produced in the 
point contact mixer and serves as a stable reference frequency 
close to the HCN laser frequency. 

The HCN laser frequency and a second klystron, K^, of any 
convenient frequency are also applied to the mixer resulting in 
the laser frequency plus and minus a sideband. The difference 
frequency between one of the sidebands and the stable reference 
is applied to a second lock box which electronically tunes K 2 . 

At this point,. K 2 exactly compensates for any drift in the laser, 
producing a stable sideband offset from the stable reference, but 
phase-locked to the frequency standard. The HCN laset and 
outputs are then applied to the next mixer and their 12th harmonics 
are used as a stable reference falling near the 28y H^O laser 
frequency. 

In a similar manner, a chain of stable sidebands may be 
constructed extending to near optical frequencies. The upper 
frequency limit to such a chain 'is given by the frequency res- 
ponse of the mixer elements. Previous work of Sokoloff, Sanchez, 
Osgood and Javan^ has demonstrated the response of the point 

3 

D. Sokoloff, A. Sanchez, R. Osgood & A. Javan 
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contact diode mixers down to wavelengths near 5y, which overlaps 
the response range of optical mixers employing bulk nonlinearities 
in solids. 

Progress to date has seen the reconstruction of the HCN and 
H 2 0 lasers in order to improve their power output and stability. 
Current effort is proceeding toward optimizing the 9 . 3y and 28y 
mixing. In the next several months, a Cesium beam atomic frequency 
standard is expected to arrive on loan from the Electronic Com- 
ponents Laboratory of the United States Army, Electronics Command. 
The laser chain will be frequency locked to the Cesium standard. 
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Block diagram of a phase locked laser frequency multiplier chain. 













V. Unidirectional Laser Amplifier with a Built-in Isolator 
As described in the last progress report, an attempt is 

1 2 

being made to utilize the laser-induced line narrowing effect ’ 
to produce a unidirectional laser amplifier in a ring laser con- 
figuration. This laser would have the unique property of attenua- 
ting any feedback signal which was sent back into it, thereby pro- 
viding isolation from its surroundings. 

The N 2 0 pumping laser has been optimized to give more than 
6 W of output. The laser stability has been checked and found to 
be satisfactory by looking at the inverted Lamb dip in NH^, and 
measurements of absorption and saturation in NH^ have been carried 
out. Detailed calculations of the gain in NH^ and other systems 
are underway to determine the most suitable substance and operating 
conditions. A paper describing the operating principles of this 
device is in preparation. 
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VI . Molecular Studies Using Standing Wave Saturation Resonances 

The study of the standing-wave saturation resonance in low 
pressure CO^ 'is well underway. The effect, first observed by Freed 
and Javan, ^ is being studied in detail to determine the influence 
of collisions and intense laser fields on the detailed line shape. 
In the experiments being performed, the low pressure C0 2 absorption 
cell is placed outside the laser cavity and the CO 2 laser, which is 
spaced with 4 invar rods for rigidity and thermal stability, is 
mounted firmly on a half-ton steel rail to minimize vibration. To 
further minimize laser jitter we have found it necessary to run the 
laser under very slow flow conditions . This has been done without 
the loss of laser intensity as would have been the case with a 
sealed off laser. • Data is' being taken on several of the C0 2 laser 
lines to determine the pressure shifts and saturation parameters of 
each transition. The signal to noise obtainable bv this technique 
for locking, a laser to this narrow resonance is demonstrated in 
Figure 2. 


^ * C. Freed and A. Javan, Appl. Phys. Letters 17, 53 (1970) . 
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Fig. 2 

Standing wave saturation resonance in C0 2 gas. Full width 
of resonance is ~ 1 MHz. 
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VII. Standing Wave Features of Laser Induced Line Narrowing Effects 


A theoretical investigation has been carried out of the effect 
of an intense standing wave laser field on a coupled transition 
when the laser is tuned to the center of its atomic gain profile. 
The calculations predict a fine structure in the response of the 
coupled transition which should have important applications in high 
resolution spectroscopy and laser frequency stabilization. This 
work has been submitted for publication and is included as Appendix 
D. 


- 16 - 


APPENDIX A 

V — 

v * 

V . 

HYPERFINE STRUCTURE OF EXCITED STATES AND QUADRUPOLE MOMENT ' 

OF Ne 21 USING LASER INDUCED LINE NARROWING TECHNIQUES* . '' 

T. V7. Ducas, M. S. Feld, L. W. Ryan, Jr. , N. Skribanowitz and A. Javan 


INTRODUCTION 


( 1 ) 


Laser Induced Line Narrov/ing Spectroscopy'’ has been used to 
observe optical hyperfine structure in Ne . This technique repre- 
sents an advance in several ways over other methods of investigating 
hyperfine and other closely-spaced structure in Doppler broadened 
systems. The best conventional optical spectrometers, even if not 
instrument-limited in resolution, are unable to resolve hyperfine 

transitions whose structure is buried within the Doppler profile. 

% 

Atomic beam studies in rare gases eliminate the Doppler effect, but 
have been generally limited to metastable levels, the only states with 
hyperfine structure that have lifetimes sufficiently long for beam 
techniques. Furthermore, ordinary beam techniques are incapable of 
measuring isotope shifts. 

In the present approach, an optical technique is used in 

which the Doppler width is effectively eliminated by laser-induced 

velocity selection. Measurements of the narrowed spectral lines have 

20 21 22 

yielded values for the isotope shifts of Ne • ' at 1.15y and 

0 21 
603PA, and hyperfine constants of Ne ’ , v.’hich lead to an accurate 

2i 

value for the Ne guadrupole moment. 
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The level scheme under study (Fig. 1)- consists of the hyperfine 

sublevels of the 2s^, 2p^ , and ls^ fine structure levels^ of Ne^ 

0 

(nuclear spin I = 3/2), which form the 1.15p (2s 2 - 2p^) and 60-96A 
(2p 4 - is ) cascade transitions. The magnetic dipole and electric 
quadrupole hyperfine interactions split each fine structure level 
(angular momentum J) into a number of hyperfine components, of energy 


K 


W p = Wj + A 2 


+ D 3/4K(K+l) - I(I+1)J(J+1) 
21 (21-1) J (2J-1) 


( 1 ) 


where K = F(F+1) - 1(1+1) - J(J+1) 

and F, the total angular momentum, can take on the values: 

F = I+J, I+J-l, . . . , 1 1- J [ . 


VL 


is the total energy of the h.f. level, Wj is the energy of the 
fine structure level; A and B are the magnetic dipole and electric 
quadrupole interaction constants of level J, respectively. 


LASER INDUCED LINE NARROWING EFFECT 

• A' comprehensive 'treatment of the laser induced line narrowing 

effect is given in Ref. 1. The following simplified discussion will 

provide a sufficient background for an understanding of the present 

experiment. Let us consider a particular coupled three-level system 

in Ne^ consisting of a pair of hyperfine transitions, one at 6096A 

and one at 1.15y, which share a common level. If a 1.15y traveling- 

21 

wave laser field is incident upon a sample of Ne " , the atoms over a 
narrow velocity range which are Doppler shifted into resonance with 
the applied field couple most strongly to it. This causes selective 
changes in the level’ populations over the narrow velocity range. Those 


changes have a Lorentzian lineshape centered about v = (ft-w 0 /K„) 

Z Lt Z Z 

and the velocity range is Av ~y/K where ft is the frequency of the 
1.15;i hyperfine transition, the corresponding propagation constant 
and v is the homogeneous linewidth of the laser transition. 

This change in the velocity distribution manifests itself in a 
change of intensity with a Lorentzian lineshape (change signal) in the 

O 

coupled 6096A fluorescence over a narrow section of the Doppler pro- 
file emitted along the axis of propagation of the laser field. The 
peak of this spontaneous emission signal in the forward direction (i.e. 
parallel to the propagation vector of the laser field it L ) occurs at 
ft,., = w^+k^v , an< 3 the backward (anti-parallel) change signal occurs 
at ft^to^ “ k-^v^ ; where is the atomic center frequency of the coupled 
6096A transition. v ' 

The locations of these change signals can be rewritten simply as: 


n P = “i + 'W-dr - ! 




= w i - • 


(2a) 


(2b) 


Thus the separation between change signals due to the pair (j,k) of 
closely spaced three-level systems is given by: 


ft F (j)-ft p (k) = [u 1 (j)-u 1 (k)] - [w 2 (jj-a) 2 (k)]j~- ; (3a) 

0) 

ft B (j)-ft B (k) = tu 1 (j)~(0 1 (k)] + [a) 2 (j)-0J 2 (k)]~ . (3b) 

21 


The' hyperfine interactions in Ne 


, which produce 18 such three- 
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T 


level systems, can be analyzed from Eqs. (3). Note that forward and 
backward patterns contain information on the h.f.s. of both 1.15y and 

o . • / 

609GA transitions. Separate observation of forward and backward sig- 

( 5 ) 

nals provides two distinct patterns which must be fit by a single 
set of parameters. These parameters, the hyperfine A and B constants 
for each of the three levels, can be extracted from the separations 
between the features of the spectrum. The fact that the relative posi- 
tions of the change signals are independent of laser tuning [Eg. (3)] 
greatly simplifies the analysis. 

The complete theoretical expression describing the laser induced 
change signals is given in Table I. • 

EXPERIMENTAL ARRANGEMENT 

The experimental set-up is shown in Fig. lb'. An intense single 
mode^ 8 '^ 1.15y He-Ne laser (see below) is locked to the resonance of 
a stable passive Fabry-Perot interferometer by standard methods. The 

•laser output beam is focused into an external sample cell containing 

21 ° 

Ne • at a low pressure (^0.1 torr) . The 6096A spontaneous emission 

from the sample cell, emitted in either the forward or backward ^ 
direction, is analyzed using a pressure-scanned Fabry-Perot interfer- 
ometer with a free spectral range of 4090 MHz and a finesse of 40. 

The entire apparatus is shock-mounted to reduce vibration effects. 

Linear scanning is achieved with a system of capillary tubes (11) bleed- 
ing dry air into a cannister containing the interferometer. A chop- 
per, placed between the laser and the sample cell, enables us to use 
a. lock-in amplifier to subtract off the spontaneous emission background 
and improve the signal to noise- ratio. 


In this experiment it is essential to have the laser operate 
on a single mode, since each extra mode would. produce an additional 

set of change signals. It is also important to maximize the power ir^ 

* . 

this mode since the sizes of the change signals are essentially pro- 
portional to the laser field intensity. This is particularly criti- 
cal in this experiment because the sample cell is outside the laser 
resonator. Neither short laser cavities nor conventional Michelson 
mode-selection techniques can satisfy the single mode power require- 
ment. V7e have employed instead the following arrangement:^'^ 

Within the laser resonator (length 1.6 m) along with the He-Ne 

amplifier cell (total pressure ^2 torr) we have placed an absorber 
cell filled with neon at a low pressure (^0.2 torr). The cells are 
of comparable length and the small signal gain of the amplifier is 
sevei'al times as large as the loss of the absorber. While the linear 
gain of the system is dominated by the amplifier cell, the absorber, 
which saturates more readily, contributes appreciably to the satura- 
tion behavior. The result is that the system is driven into the 
strong coupling regime, where the modes compete heavily with one another 
Thus, laser oscillation at one mode tends to suppress oscillation at 
other modes. The net result is a continuously tunable single mode 
output with an appreciable fraction of the full power of multimode opera 

tion. In our experiment we. were able to convert over 50% of the multi- 

* " 

mode power ( modes) into a single mode. 
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HYPERFXNE THEORY 


The electronic configurations of excited states in neon consist-' 

of a 2p hole coupled to an excited electron surrounding closed shells. 

_The wave functions of the resulting fine structure energy levels, ob- 

( 12 ) 

tained from the f.s. energy splittings , may be expressed in terms 
of the admixture of the appropriate LS wave functions ^ 3 ' ^ ^ . The 
A and B interaction constants which determine the hyperfine spectrum 
[Eq. (1)] may.be obtained from these wave functions, the nuclear di- 
pole moment y, and the nuclear quadrupole moment The. resulting 

expressions for the A's and B's of the three levels of interest con- 
tain the coefficients of the LS wave functions and the parameters: y; 

Q; ~ r for the 2p hole; — for the 3p electron; and the contact terms 

r r 3 

A-, and A. for the 3s and 4s states of the excited electron. 

3s 4s 


In a previous experiment 


(16) 


Grosof, Buck, Lichten and Rabi 

21 


measured the A and B constants of the Is level of Ne to high (0.1%) 

• J 


accuracy. But since 


can only be estimated to 10% accuracy, 


2p hole 


using the Fernti-Segre formula 


(15) 


and their measured value of A(ls c ), 

D 


their value of Q has this same uncertainty [B(ls c ) Q~. — 

d r 


We make use of the measured 

(17) 


(16) 


2p hole 




values of A (Is..) and B(ls c ) and 

D D 


of the known'*' ' value of y to reduce the number of parameters in 

our theoretical description by three: Q is expressed in terms of 

B (] s ) and ~ 3 - ; and A- may be expressed in terms of A(ls_) 

5 r 2p hole ... 

i-> (independent of the admixture of states for the ls c level). 

r 2p hole 5 

1 . 1 


The remaining parameters are: — 3 

r 2p hole 


and A 


4s 


3p electron 

The contribution of the 3p electron to the h.f.s. is negligible. Fur- 
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thermore, the error in our estimate of (from the Fermi-Segre form 

ula) produces only a small uncertainty in A(2s 2 ). Thus, the h.f. - 

interaction constants are sensitive only to the parameter 

2p hole 

In the following paragraphs we expand upon the preceding discussion 
and describe the particular calculations involved in obtaining exprcs 
sions for the hyperfine constants. 

In our analysis we use the wave functions of Vainshtein and 

M- (14) 

Minaeva : 


l? s 2 “ 

-.660 | 

I 3 ?! + 

.752] 

I 1 ?! / 


1 2p 4 *= 

.127 

\\ - 

. 309 

| 1 D 2 + . 9 4 3 | 3 p 2 , 

(4) 

|1® 4 = 

.964 

l 3 Pl ♦ 

.266 

,1 

1 Py • 



Where 2s+1 L T are wave functions in the LS representation. These 

J 

wave functions have been derived from the fine structure energy 

splittings ^ 2 ^ , taking into. account electrostatic, spin-orbit, spin- 

spin and spin-other-orbit interactions. 

. The LS matrix elements for the magnetic dipole and electric 

5 5 1 

quadrupole interactions for the p 3 and p p configurations are 
given by Childs ^ ^ . We use the non-relativistic limit for these 
matrix elements since relativistic corrections are negligible (<1% 
for low Z atoms such as neon) combining these with Eqs. (4) we get 
the following expressions for the A and- B interaction constants for 
the three levels of concern: 


A(2s 2 ) 

= 1.352 

a 2 P 

hole + 

.460 

a. 

4s, 

B(2s 2 ) 

= -.139 

b 2p 

hole' 



A(2p 4 ) 

= ..367 

a 2p 

hole + 

.403 

3p 

B(2p 4 ) 

= .053 

b 2p 

hole 

.219 

b. ; 
3p 


A(1s 4 ) .819 a 2p hole + .051 a 3g ^ 

B(1s 4 ) = .158 b 2p hole ? 


With 


_ 2puB 1 

np I r 

. 2Q 1 

b = e — 

. n P ' r3 np' 


np 


where y D is the Bohr magneton. 

Nov/, if we employ the results of Grosof et al 
level of Ne 2 ^: 


(16) 


for the Is, 


A (ls c ) = -267 . 68± . 03MHz 

D 

B(1s 5 ) = -111 . 55± . 10MHz 


(18) 


and the fact that 
wave function: 

. A (Is,.) = 1/4 a_ + 2/5 a 


, independent of the LS coefficients of the Is 


'5' ~ / ‘ 3s 

B(ls_) = -2/5 b. 


2p hole 

5' ~ / ” ~2p hole 

v/e can express our A and B constants in MHz in the following manner 


A(2s 2 ) = .460 a 4s - 843.27 R h , 

B ( 2 s 2 ) = -38.786; 


A(2p 4 ) = -228.78R h - 251.39R e , 
B(2p 4 ) = 14.898-61.191 R e /\; 


( 6 ) 


where 


A(ls 4 ) = -54.639 - 459.74R h , 
B(1 s 4 ) 43.940; 


2p hole 


= R^ x 10 26 cm 3 , 


pr = R x 10 26 cm -3 , 

,3p electron 


Equation (6) provides quantitative expressions for the con- 
tributions of the different parameters to the A's and B's. Esti- 
mates of A. and R can be obtained from the fine structure energy 
4s e 

(15 ) 

levels using formulae in Kopfermann . The uncertainties in these 
calculations are of the order of 10%, but the uncertainties they intro 
duce into A(2s 2 ), A(2p 4 ) and B(2p 4 ) are much smaller, on the order of 
1%. This is so becuase R e is less than 1% as large as R^ and, as 
will be seen below, the A^ s term in Eq. (6) contributes less than 10% 
to A(2s 2 ) . ' . . 

The preceding discussion indicates that the single dominant in- 
fluence in determining each hyperfine interaction constant is 

i-r . Because of the sensitivity of our expressions to this 

r 2p hole 

parameter, we must provide for its variation with the state of the 
excited electron. Equations (5) and (6) can be modified in a straight 
forward way to include such effects. The B of each level is less 
than 10% of the corresponding A in magnitude, so that small changes 
in 1/r^ 2p h 0 ]L e ^ ave negligible influence on the h.f.s. through the 
B's. (This assertion is born out in our computer fits to the data.) 
But these modifications cannot be ignored in the A's. The A's from 
Eq. (6) are therefore written: 

A(2s 2 ) = .460 a 4c . - 843.27R h (2s 2 ), 

A(2p 4 ) = -228.78R h (2p 4 ) - 251.39R e , * (7) 

A(1s 4 ) = -54.639 - 459.74R h (ls 4 ) . 
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DATA ANALYSIS AND THEORETICAL FIT 

Figure shov/s typical experimental spectra taken with the 
Fabry-Perot interferometer, and the computer-generated fit for for- 
ward (Fig. 3a) and backward (Fig. 3b) runs. In each case the upper 
experimental trace shov/s the narrow change signals, and the lower, 
trace shows the Doppler-broadened background observed .without the 
lock-in amplifier. Note that the hyperfine structure is ordinarily 
completely masked by the Doppler widths. 

The laser frequency, which was carefully monitored, was stable . 
to about 2MHz over the course of a run (about one minute) . This 
figure is considerably smaller than the linewidths of the observed 
change signals, which are of the order of 250 MHz (including a Fabry- 
perot width of 'VlOO MHz) . The largest change signals in Fig. 3 have 
intensities of about 0.5% of the intensity of. the spontaneous emission 
background. It must be emphasized that the observed spectra are 

O 

richer than ordinary 609 6A h.f. spectra in that our data contains in- 
formation about both -the laser transition and .the coupled spontaneous 
emission transition. The markings on Fig. 3a denote the eight 3-level 
systems producing the observed change signals (the contribution of the 
remaining ten change signals being negligible) . The two additional 

features in the forward and backward traces are due to the presence 

?0 22 2 ] ( 20 ) 
of Ne and Ne in the sample of Ne 

A PDP-12 computer was used to display intensity patterns employ- 
ing the expressions in Eq. (1) and (2) and in Table I. Trial values. 

for the A 1 s and B's were obtained from Eq. (6), using Fermi-Segre esti 

3 3 

mates for A. and 1/r - and an optimum value of 1/r „ . , . Only 

4 S 3p ' 2p hole 

small deviations from these trial values were required to find the 


TT 


hyperfine constants that provided the best fit to the data. These 

small deviations from the values generated from Eq. (6) can be 

associated with the theoretical considerations which led to Eq. (7), 

• *. 

3 

where a slightly different 1/r ^ ho] e characterizes each A. 

The *a(j)'s and B(j)'s appearing in the intensity formula (Table I), 

which determine the relative intensities of the change signals, were 

(13) 

calculated from the electric dipolematrix elements . A convenient 

simplification is that the a(j)'s and 3 ( j ) 1 s are independent of the 

\ 

admixture of the LS wave functions for a particular fine structure 
energy level (aside from a proportionality constant) . 

Sixty-six runs were analyzed in all: 41 forward and 25 back- 

ward. In analyzing the data we operated under the strict require- 
ment, imposed by the theory of laser induced line narrowing, that the 
same set of A and B values has to fit both forward and backward ex- 
perimental traces. 

In several experimental runs, the laser frequency was intentionally 
shifted vlOOMHz. This alters the intensity pattern in prescribed 
ways (Table I) , serving as a further check on our identification of 
the features. 

The best fit to the data (Table II) provides values of : 

R^(ls^) = . 865 ± .015; (2p^, ) = . 948; R h^ s 2^ = .866. The slightly 

larger error in the latter value is due to the uncertainty in esti- 
mating a^ g . 

3 . 

Combining 1/r hole ^ or ^ s 4 ^ eve ^ w ith B(ls^) from Ref. 18, 

we obtain an accurate value for the uncorrected nuclear quadrupole 

21 

moment Q' for the ls^ level of Me . This is in gopd agreement with 
the value quoted in Ref. 18. Estimates of Q' for the other two levels 


C_ 


would be less accurate since our experiment is not very sensitive to 
the precise values of the corresponding B's. To establish the actual 
Q, the Sternheimer correctons would have to be taken into account. 

We can also obtain a value for a^ s » t * le Ferm i contact interaction 

3 

„ constant, from Eq. (5) and our value of 1/r 2 ^ hole ^ or t ‘ le ^ S 4 level 

Note that this value if 15% smaller in magnitude than that calculated 

from the Fermi-Segre formula. 

3 

Our values of 1/r 2p hole ^ or ^ S 2 an< ^ ^ S 4 ^ eve -*- s ' which 

have the excited electron in an s-state, are within 1% of each other. 
The corresponding result for the 2p^ level is about 10% larger than 
these. The former two values are in close agreement with the semi- 
empirical value for the Is,, level from Ref. 18. 

It is interesting to compare these values with that for the neon 
ion ground state, in which the excited electron surrounding the 2p 
'hole is removed. A calculation using restricted Hartree-Fock wave 
functions gives a value • for 1/r 3 2p ho] e which is within 2% of 
our results for the 2s 2 and ls^ levels. 

. These conclusions are, of course, subject to the approximations 
made in deriving the A's of Eq. (7) and the B's of Eq. (6). One 
source of error is the uncertainty in the LS admixture coefficients 
of our wave functions [Eq. (4)]. These values .are expected to be 
accurate for the ls 4 and 2s 2 levels, but somewhat less certain for 
the 2p 4 level. This is primarily due to the difficulties associated 
with properly including effects from spin-spin and spin-other-orbit 
interactions^ 22 ^. Note that the value of A(2p^) is particularly 
sensitive to small deviations in the admixture coefficients, since 
the 2 d^ state is ^89% a pure LS state. 
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The assumption of a single 1 /r ^ p hole c h aract erizing each 

(23) 

level ignores core polarization effects ' , which result from the 

2 2 

spin-dependent distortion of the (Is) ( 2 s) core by the outer elec- 
trons. Taking this effect into account theoretically would necessi- 

3 

tate for each level: the use of different values of 1 /r hole 

in the magnetic h.f.s. resulting from the orbital and spin moments of 
the electrons; a separate value of 1 /r 2 p hole assoc i ate< ^ with the 
quadrupole interaction; and would also add a small contact term to 
the expressions for the A's. 

It should be possible to explore fully the details of the above 

considerations by studying the laser induced change signals produced 

at other transitions branching from the upper and lower levels of 

the 1.15y and other neon laser lines. 

209122 

We have also measured the Ne 11 isotope shifts at 1.15y and 

° 20 2 2 * 

6096A. Ne and Ne were introduced into the sample cell containing 

on , 

Ne to bring up the intensity of their change signals, thus facili- 
tating their identification and measurement. The shifts for the 1.15y 

O 

and 6096A transitions are extracted from the change signal spectra 

21 

using Eq. (3) [and Eq. (1) in the case of Ne ], and are summarized 

20 22 

in Table III. The Ne ' ‘ results are in excellent agreement with 

/ 2 04 2 5) 

previously measured values' ' ' . Note in particular that the 1.15y 

transition of Ne is shifted to a higher frequency than that of 
22 

Ne . An analysis of the results is in progress and will be reported 


later. 
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CONCLUSION 

* . 

We have demonstrated the applicability of laser induced line 
narrowing techniques in observing optical h.f.s. in a system whose 
features are ordinarily masked by Doppler broadening. These tech- 
niques can be applied in principle to measurements of closely-spaced 
structure of transitions coupled directly or indirectly through 
radiative cascade to either of the laser levels. 

Although the change signals in this experiment are considerably 
narrower than the Doppler width, for many other systems they may be 
two or three orders of magnitude narrower still. The limits on the 
resolution are only determined by the homogeneous width of the lines 
and the resolution of the detection apparatus. 
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where E refers to the sum over all possible cascade 

j 

transitions of the type F-’-F’-dT 11 

|ct ( j ) J 2 = |m x (F , # M f 1 jF 1 , ,M f m ) | 2 + |y z (F 1 ,M F 1 ;F 11 ,M F 11 ) | 2 , 

O 

• • for observation of both polarizations of the 6096A 
radiation together. 


| 3 ( j ) | 2 = |v_(f’,M f ;F 1 ,M p 1 ) I 2 ; 


p's are the hyperfine electric dipole matrix elements and 
can be obtained from tables in Ref. 13. 

Note that the relative magnitudes of the a(j)'s and 3 ( j ) 1 s 
are independent of the fine structure coupling of the levels 

O 

is the frequency of the spontaneous emission at 6096A. 

ftp are defined in the text 
B 

Explicit expressions for the proportionality constants 
and C 2 are given in Ref. 1. 

The line-widths of the change signals (full width at half- 
maximum) are: 

\ = *1 + >4 <V Y 2 ,±Y 0> 

where y^ is the decay rate of level i. For our case i=0,l,2 
correspond tb the 2p^, ls^ and 2s 2 levels respectively. 

C 

Sec Ref. 1 for details. 


TABLE I 
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Ne f.s. levels-*- 

1S 4 ‘ 

2 p 4 

2s 2 





Magnetic Dipole h.f. 
Interaction constant, 
A (MHz) 

-452 

-.219 

-756 

Electric Quadrupole 
h.f. interaction 
constant, B (MHz) 

+ 44 

+ 15 

- 39 

R h ( 1/r3 2p hole 
in 10 2 6 cm" 3 ) 

. 865+. 015 

t • * 

.948 

.866 


Nuclear Quadrupole Moment : 

Q' (ls 4 ) = .0926+. 0016 b* Q'(ls 5 ) = .093+.010 b 

Fermi Contact term f o r 3s electro n : 

(theoretical) 

a_ = -207 MHz * . a- = -245 MHz *** 

3s 3s 

2p hole 

R h (ls 4 ) = . 865+ . 015* R h( ls 5 ) = * 86 ± 

R h (Ne + ) = .84‘ 

- * Present work 

** Ref. 16 • 

*** Fermi -Seg re Estimate 

Semi-Empirical value derived from Ref. 16 
Ref. 21 


TABLE II 


v(Ne 22 ) - v (Ne 20 ) 


• 

'I 

Present 

Work 

Ref. 

21 

Ref . 
23 

Ref. 

24 

o 





6096A’ 

1708 

1706+30 

1680 + 


1 . 15y 

261 

257+8 


261+3 


v (Ne 22 ) - v (Ne 21 ) 


6096A 


1. 15u 


Present Work 


1010 


-186 
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FIGURE CAPTIONS 
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Fig. 1 Hyperfine splitting of ls^, 2p 4 2s 2 levels in Ne 

Fig. 2 Experimental Arrangement. Note the different configurations 
for observing forward and backward 6096A spontaneous emis- 
sion. 

Fig. 3 Comparison of experimental traces and theoretical fit. Note 
that the Doppler profiles for forward and backward traces 
completely mask the hyperfine structure. The change sig- 
nals labeled are the eight most prominent, and are associa- 
ted with three level cascade systems of the type F->-F 1 ->F 11 . 
Denoting these systmes as ordered triplets, they are: A= 

(3/2, 3/2, 1/2); B= (3/2, 5/2, 5/2); C= (3/2, 3/2, 3/2); 

• D= (5/2, 7/2, 5/2); E= (5/2, 5/2, 3/2); F= (1/2, 1/2, 1/2); 

G= (3/2, 5/2, 3/2); H= (5/2, 5/2, 5/2). 

For these computer fits to the data y =265 MHz and y=225MHz. 

r 

Of these widths ^lOOMHz is instrumental. Also, the Doppler 
width k^u=2500MHz. 

TABLE CAPTIONS 

Table 1 Theoretical expression for intensity pattern of change signals. 

Table 2 Experimental .results for h.f. interaction constants and 

values for 1/r* , , , O', and a_ . 

' 2p hole 3s ' 

Table 3 Observed isotope shifts at 6096A and 1.15y in Ne^*^'^'* ' 
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APPENDIX B. OBSERVATION OF INTENSE SUPERRADIANT 

EMISSION IN THE HIGH GAIN INFRARED 
TRANSITIONS OF HF AND DF MOLECULES 


J. Goldhar, R. M. Osgood, Jr. and A. Javan. 


ABSTRACT 

Intense stimulated emission (superradiance) is 
obtained in several high gain transitions of HF 
and DF molecules. The narrowing of the linewidths 
of the high gain transitions are observed in detail. 
Pulsed energy and peak powers are comparable with 
those obtainable in a similar transverse CO^ laser. 


2 . 


This letter reports the observations of intense stimulated 
emission (superradiance lin a number of high gain, rotational- 
vibrational transitions of HP and DP molecules in the 2.7y and 
3.8y range of wave-length respectively. The vibrationally 

excited molecules are produced with a transverse pulsed discharge 
in a flowing mixture of molecular hydrogen (or deuterium) and 
SF, gas, at relatively high pressure. The spectral distribution 
of the individual high gain transitions have been studied v/ith a 
pneumatically tuned Fabry-perot interferometer. Each of the 
superradiant transitions is found to show considerable line-width 
narrowing due to substantial amplification of the traveling 
optical wave. In addition, it is found that in the presence of 
a regenerative optical feedback (i.e . with the device used as a 
laser with an optical resonator) , the laser oscillation on each 
line occurs mainly in a single resonator mode. With regard to 
the obtainable peak power and the pulse energy , the performance 
of the system is competitive with that obtained at 10. 6y using 
a CO 0 laser excited by means of a similar transverse pulsed 
discharge. 

In the experiment, the laser configuration consisted of a 
50 cm long transverse discharge with 101 equally spaced pm 
electrodes arranged in a row facing a long cylindrical brass, 
anode across a 2.5 cm gap. A 0.01 yF capacitor at high voltage 
is discharged across the tube through a 1000 ohm resistor in 
scries with each pin. The discharge tube was terminated with Ca ?2 
Brewster windows. Gases were introduced into the system by 
means of standard needle valves and flow 7 was maintained with a 


10 CPM fere pump. 


In the experiment, most of the detailed observations were 

made with the hydrogen gas. Typically, the partial pressure 

ratio of SF^ to was 10:1. In this system, the vibrationally 

. 4 . 

excited H7 molecules are formed through chemical reaction of 

the fluorine atoms, (which exist as a byproduct of the SF g dis- 

2 

sociation) , and the molecular hydrogen. As noted earlier, 

this reaction is 31.7 kcal/mole exothermic for HF and DF, and 

hence capable of producing vibrationally excited HF up to the 
vibrational level 
v=3/(and up to v=4 in DF) . 

The distribution of infrared output radiation among tne 
various rotation-vibration transitions in HF , is found to be 
dependent on the operating conditions such as flow rate, partial 
and total pressure, and applied voltage. If the device is 
operated as a laser with a regenerative optical feedback, at 
relatively low pressures of SF^ and ^ and low applied voltage, 
(below about ICk.v. ), the output is distributed among several 
rotation-vibration transitions of each of the three v=3->2, 2-*l 
and 1^0 vibrational bands. ^ At the elevated pressures and ror 
the applied voltages in the regions of about 15 kilovolts (or 
higher) , the laser output spectrum tends to peak in the v=l->0 
transition with P (4) the most energetic line. This behaviour 
indicates that at the elevated pressures and voltages, (where the 
HF population is appreciable) , the HF molecules in the v-*l states 

decay rapidly to the v=l state, causing a build-up or population 

, 6 

in the v-i s cate , which apparently decays at a slower rare. 

An important characteristic of this system is its capability 
of producing very high gains leading to intense stimulated 


scission without regenerative optical feedback (superradiance) . 

In analyzing the output spectrum of the superradiant device, it 
was most essential to guard against unwanted optical feedback 
introduced by small reflections or scattering from auxiliary 
optical components present in the path of the output beam. This 
was done by introducing heavily absorbing filters placed at an 
angle in the path of the output beam and appropriately mis- 
aligning the various components to prevent the feedback. Such 
precautions were particularly important when a folding mirror 
was placed at one end of the amplifier to obtain double pass 
amplification at the output. Fig. 1 shows the method used to 
monitor the presence of an unwanted optical feedback. A sample 
outpuc was provided by placing a beam splitter on the amplifier 
axis near the folding mirror. When the light emanating from 
the output end of the amplifier was interrupted, the presence 
of any feedback exterior to the amplifier could then be detected 
as a signal appearing at the sampling output. 

In the high pressure region (about 200 torr SFg and about 
20 torr H 0 ) , it was possible to obtain intense stimulated emission 
on mhe P(4), v=l-*0 and P(3), v=2+l transitions with the amplifier 
used in a single pass; the P(4), v=l-»-0 was generally found to be 
more intense. 

A pneumatically scanned Fabry-Perot interferometer with a 
S.25 cm mirror separation was used to analyze the spectral dis- 
tribution of each high gain transition. In this system, the 
ouuput of the interferometer was detected through a small pinhole 
wiah an infrared detector. The interferometric analysis of the 
sooctral distribution in each superradiant transition showed a 


5. 


ruxj. 


W CL S 


This 


-width at half maximum below 110 MHz. The 110 MHz limit 
instrumental, determined by the finesse of the Fabry-Perot. 
result was obtained with the amplifier used in both the 


single and double pass configurations. 

An additional noteworthy behaviour of the system is that 

under the superradiant conditions, it is possible to obtain 

stimulated emission with sufficient intensity to cause amplifier 

saturation even for a single pass amplifier. This was observed 

by noting that at the limit of highest gain, the double-pass 

transition 

amplifier output on the P(4), v=l-»0^was only about 20 times 
larger than that obtained in the single pass. Furthermore, using 
the feedback monitor described above, the sampled signal of the 
double pass amplifier with no feedback was compared with that of 
uhe amplifier used as a laser with a complete regenerative feed- 
back. The power level in the latter case was found to be only a 
few times higher than that of the double pass amplifier 

Under the high pressure condition where the amplifier gain 
is near optimum, the spontaneous emission linewidth is essentially 
due to collision broadening. The exact magnitude of this 
broadening is not as yet known for the HF transitions. However, 
in the presence of about 200 torr SF„ pressure, the spontaneous 

b 

emission line width car. be assumed to be at least about 3000 MHz. 

observed corresponds to 

Accordingly, the* 110 MHz line width limit /, a line narrowing 

of 30 or greater. In the absence of saturation, this factor is 
known to be given by the square root of GL v/here G is magnitude 
of the gain per unit length and L is the amplification path 
length. However, it is important to guard against hasty appli- 
cation of this relationship in obtaining an accurate measure 


6 . 


of the amplifier gain in this system. For instance., it is 
probable that the presence of each pin electrode in the amplifying 
medium can cause appreciable scattering of the infrared emission 

to introduce sufficient feedback to cause additional narrowing 

7 

of the line profile. Furthermore, the saturation effect is also 
expected to play an appreciable effect on the line-shape. 

The frequency spectrum of each oscillating line was also 
analyzed when the device was used as a laser with two aligned 
mirrors providing complete regenerative feedback. Interestingly, 
the frequency spectrum of each oscillating line, as averaged 
over a number of pulses, was found to be distributed within a 
frequency interval below about 70 MHz resolution limit of the 

g 

Fabry-Perot . In this case, a much broader spectrum is expected 

because of the possibility of the laser oscillating on several 

longitudinal modes of the laser resonator which was 11 8 cm long 

(corresponding to 130 MHz mode frequency spacing) . The narrow 

line-width for each oscillating line was observed while the gas 

pressure and the applied voltage and hence the gain in the 

various transitions were varied over a wide range. (The lowest 

somewhat 

total gas pressure was /less than 10 torr) . While it is tempting 
no attribute this persistent behaviour to the high gain property 
of the medium, it is more likely that the effect originates from 
mode coupling in the presence of collision broadening. For a 
collision broadened line, the saturation effect decreases the 
gain uniformly over the whole line profile leading tc mode 
competition . In general, however, complete single mod in g is 
not always expected for a collision broadened line due to the 
so called ” spatial hole burning" effect. This effect occurs 
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in a region of space within the resonator where the crest of the 
standing wave for one mode falls near the node of another. This 
would allow coupling of both modes to the amplifying medium when 
the molecular mean free path is appreciably below one wavelength 
- a condition valid for a fully collision broadened line. 

In our system, however, the spatial hole burning of this 

type may be smeared because of rapid density fluctuations and 

the acoustical waves caused by the high current discharge pulse . 

This would inhibit multimoding on, say, several longitudinal 

modes but would allow simultaneous oscillation on a longitudinal 

and a transverse mode having appreciably different intensity 

distribution in the direction transverse to the laser axis. 

The latter possibility was verified by using an appropriate pair 

of curved laser mirrors to introduce transverse modes separated 

in frequency by tens of MHz from the corresponding longitudinal 

modes. With a fast response infrared detector, it was possible 

to obtain beat notes in the appropriate frequency regions 

dependent, as expected, on the degree of aperturing of the laser 

simultaneous 

resonator. This result indicates the possible^/presence of one 
longitudinal and several transverse modes at frequencies falling 
within the 70 X Hz limit of the Fabry-Perot. 

Study of the output power obtainable from this sytem has 
revealed an important practical property; namely, that the system 
is potentially capable of supplying considerable energies and 
high peak powers in each pulse. The optimum power and efficiency 
was found to be an increasing function of the applied voltage. 

At the 20 kv limit of the available high voltage supply, it was 
possible to obtain 16 kw peak power (corresponding to 0.6% 
efficiency) in a pulse of 700 ns duration. This was obtained in 
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a mixture of S?,., K- v;ith partial pressures of 300 and 10 torr 

O ^ 

respectively. This result suggests that by increasing the 
applied voltage and scaling the volume and the gas pressure, it 
is possible to obtain much larger output powers, competitive with 
those reportedly obtainable in the transverse discharge CO^ laser. 

With the gas used instead of K 0 , it was possible to obtain 
intense stimulated emissions on the DF transitions in a manner 
similar to that observed in HF, as described above. With the 
amplifier used as a laser with regenerative feedback, the DF 
oscillations appeared over several transitions in the v=4~ > 3; 
v-3->2 , v=2->i, and v=l->-0 bancs. In the high pressure region, 
the lines capable of intense superradiant emission were P(5) 
and ?(S) of both the v=2-»-l and V=3-*v=2 bands. More detailed 
studies of the pressure dependence of the DF spectrum are now 


in progress. 

In conclusion, we would like to point out that the high 
pressure transverse gas discharge provides a general method of 
extracting high peak powers and pulse energy for many other gas 
laser transitions. For instance, we have recently observed intense 
laser oscillations in the 5y bands of Pure CO excited by a trans- 
verse gas discharge at relatively high pressures. Detailed study 
of the system will be published separately 

Wo would like to acknowledge Dr. ?. Zernike for making avail- 
able the Fabry-Perot flats used in the experiment. We would also 
wish to thank Mr. Ronald McNair and W. L. Ryan for their expert 
technical assistance. Lastly, we are grateful to Prof. S. V. 


George anc. Dr. R. Carbone for numerous helpful discussions. 
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possibility of obtaining large gains in an electrically 
pulsed, chemical HF laser. 

2. For chemical HF, DF lasers produced in a supersonic 
flowing mixture of D 0 or with thermally heated SF^, 

see D. J. Spencer, H. Mirels, and T. A. Jacobs, Appl. Phys. 
•Letters 16 , 284, 1970. For the application of flash 
photolysis, see J. K. Parkes and G. C. Pimentel, J. of 
Chem Phys., 5JL, 91, 1969. 

3. A. J. Beaulieu, Appl. Phys. Letters, 16_, 504, 1970. 

4. Parkes and Pimentel, J. of Chem. Phys., 51 , 91, 1969. 

5. The transitions in the v=3‘ > 2 and v = 2~*1 bands ranged from 
P(3) to P(6) with P(3) generally being more intense. In 
the v^l^O band the transitions P(4) through P(6) were ob- 
served with P (4) being the most intense. 

6. In. this process, the rapid decay of the HF molecules in the v>] 
states to the v=l will be dominated by v-v type collisions 
with HF molecules in the v=0 state. However, the slower 
decay of the v=l state will result from collisions with 

^*2 an< -- ^*(3 molecules. 
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-he nigh Cc.in conditions , it was possible to tricger 
.'.plifier into oscillation" with a diffuse regenerative 
:x. For instance , by holding a sheet of paper at the 
output end of the folded amplifier, it was possible to 
'dot air. strong oscillations. 

“Tnis could ce detected by means of the feedback monitor (see 
Fig. 1) . For a discussion of non-resonant feedback and its 
implications see R. V. Ambartsumyan, N. G. Basov, P. G. Kryukov 
ana V. C. Letokhov JET? 24, 481 1857; and R. V. Anbarsumvan, 

P. G. Kryukov and V. C. Letokhov JETP 24, 1129, 1967. 

'.■..isn me device was used as a laser with re Generative feed- 
back, the interferometer could be better alignfed with respect 
to me incident beam without introducing feedback into the 
laser. Accordingly, a better system finesse was then obtained. 
Ar. additional source of the smearing of the spatial hole 
burning may also arise from a rapid variation of the refrac- 
tive index during the build up of radiation within the 
resonator causing a wavelength sweeping of each mode. Such 
an index change may arise from molecular dissociation due to 


n.- 


the current pulse and 
v ..in amplifying trans 
index. 

Csgocd, Goidhar and 


the build-up of population in the high 
it ions contributing to the refractive 

McNair to be published in Journal of 


Quantum Electronics. 



CAPTIONS 


Fig. 1: Experimental apparatus for Analyzing the Behaviour 

of the Superradiant Device. 

A) Folding Mirror; B) Sampling; 

V- /i { -»v o/j .l«f ■» J 

C ) Sampling detector with filter beamsplitter; 

D) Discharge tube; E) Absorbing Filters 

F) Monochrometer s G) Collimator 

H)_ Febry-Perot Interferomater with_ pressure scanning 
valve, pinhole and detector. 


APPENDIX C. 


FLUORESCENCE INDUCED BY COHERENT 


OPTICAL PULSES 


H. P. Grieneisen, N. A. Kurnit and A. Szoke 


ABSTRACT 


A coherent electromagnetic field, resonant with 
an optical transition, drives the molecules of 
a gas into the upper and lower states alternately, 
unless relaxation is too fast. The observation of 
this coherent excitation of optical levels by fluor- 
escence measurement is discussed and experimental 
results are presented. The 4.3y fluorescence of 
C0 2 gas excited by a 10. 6y laser is studied; values 
for the dipole moment and relaxation times are ob- 
tained in good agreement with those obtained by other 


methods . 
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1 2 

Studies of photon echoes, adiabatic rapid passage, tran- 

*5 4 

sient nutation, and self-induced transparency have extended into 


the optical region many of the methods used in magnetic resonance 
for probing relaxation processes. This letter presents results 
on a related, sensitive method of examining relaxation effects 


^11J_ WillCi V_ JL ^ *** u a; 

case, a simple and direct correspondence between the analysis 
and the experimental results emerges; complications from non- 


We consider an atom or molecule in a gas with three energy 
levels of primary interest. The system is initially in state 1; 
a resonant optical pulse excites the system coherently to state 
2, and the system subsequently decays to state 3 (possibly 
identical to state 1) while radiating spontaneously. The optical 
electric field of the running' wave pulse, at the position of the 
atom, x, is described by; 


where e is the polarization vector , £(x,t) is a slowly varying 
envelope, and w»w 12 is the optical carrier frequency, possibly 
also varying slowly in time. During the pulse the equation of 
motion of the density matrix p can thus be approximated by that 
of a two- level system: 



3 4 

linear propagation effects ' ' do not appear. 


E (x, t) = e£(x,t)cos (wt-kx) 


(1) 


P = -^[H, pi- 



rn 


where the Hamiltonian H = H -y *E(x,t) consists of the unper- 


turbed part, H q , having the energy levels 1, 2, and the inter- 
action of the optical field with the electric dipole moment of 
the atom, y . The relaxation is treated phenomenologically 
(Bloch equations) , it is characterized by the decay time of the 
energy, T^, and that of the coherent dipole moment, T^. The 

g 

relaxation can be caused by radiative processes, collisions, or 

diffusion. A more detailed account is beyond the scope of this 

letter. Eq. (2) can be integrated (as a rule numerically) , with 

the initial condition j (t=0 ) = 5 j]/ to give the value 

of p. . (t=t ) after passage of the pulse. If for simplicity 
ID a 

we assume that the lifetime of the excited population is much 
longer than T 2 , the decay of P22^ t>t: a +T 2^ can be considered sep- 
arately since Thus the fluorescence intensity probes 

the upper state population at the end of the pulse. We can 
illustrate the expected results in some simple cases: 

1) No relaxation, field on. resonance for all the atoms, 

Eq. (2) can be integrated immediately in the rotating wave 
approximation to yield: 


P 22 (t a } 


P 11 Sin 


i 


(y£(x,t)/2ft)dt 


(3) 


where y = | ( e • V* 0 p ^ 1 2 ^ a dramatic effect. Increasing the 

field (or lengthening the pulse) first increases, then decreases 
the observed intensity. The period of the modulation yields the 
value of the dipole moment y. Also it is a direct measure of the 

/ +oo 

£(x,t)dt, introduced by McCall and Hahn, 4 

-CO 

which has importance in propagation effects. As an optically thin 
sample changes the value of 9 only slightly, this experiment can 
be used as a diagnostic tool in self-induced transparency and 


related studies. The angle 0 is also related to the pulse 

2 2 2 . 

energy/unit area, S, by 0 = 8ny ST^/cft , provided we define 

the pulse width as t = [ f £(x,t)dt ] 2 / / £ 2 (x,t)dt. 

P J-co ‘'-oo 

2) For fast relaxation, T i' T 2 <<T p' t * ie stea< ty state 
solution applies. For a square resonant optical pulse one has 


p 22 (t a J 


p° x (y£/h) 2 T 1 T, 


2 l+(y£/h) 2 T 1 T 2 


(4) 


It can be seen that the fluorescence increases monotonically with 
pulse energy, first linearly , then more slowly. The saturation 
field measures the product T..T.-. 

In our experiment we deal with a spatially degenerate 
rotational-vibrational transition of CO^. The measured fluorescence 
intensity is then proportional to 


p 22 (t a* 




d (Aw) g (Aco) 



+J 

.j P 2M,2M (A “'V 


(5) 


where g(Aw) is a spectral distribution function (a Gaussian of 
width AiOp for Doppler-broadening) . 9 In the case of spatial 
degeneracy, the selection rule AM=0 for a linearly polarized 
field allows the transition to be treated as a collection of two- 
level systems, each one with its own resonance frequency and 
matrix element y M (J) = y Q (J) (1-M 2 /J 2 ) 1/2 for a P-branch transition^ 0 
Eq. (2) and the integral (5) have been evaluated for width Aw^/2tt 
= 60MHz and for various pulse shapes on a digital computer 
and results are presented in Fig. 1 for some cases of interest. 

The fluorescence intensity is plotted against laser field strength 
for various T^Tj values. The salient features are: At long 

relaxation times (low pressures) there is modulation even for 


high J transitions,^ the modulation disappears at T 2 ~t^. 

The fluorescence intensity is always linearly proportional to 

the laser intensity (quadratic in £ ) at low laser powers. For 

T^,T 2 <<t it becomes proportional to the square root of laser 

l /o 

intensity (linear in £) when (y£/fi) (T^T^) 7/ >i and finally 
saturates when (y£/ft) (T 1 /T. ) ^ The parameters y,T-, 

-L Z V/ U Z 

1/2 

(TiT 2 ) can be estimated from the first minimum, the disap- 
pearance of the modulation, and the linear to square root 
transition region respectively. 

Our experimental arrangement is shown in Fig. 2a. A 4m 
flowing C^-^-He laser (A), operating on the P(20) 10°0-00°1 
10. 6y transition, with line selection provided by a diffraction 
grating (B) , is Q-switched by a rotating mirror (C) . Apertures 
(D) placed at both ends of the cavity each introduce ^30% 
diffraction loss in order to eliminate off-axis modes and improve 
line selection, resulting in a pulse with good amplitude and 
shape stability (Fig. 2b) . Frequency drift of the laser was 
corrected by moving the 90% reflecting laser output mirror (E) 
with a piezoelectric transducer (F) . The laser gain was kept 
sufficiently low that detuning of the cavity mode by more than 
MMHz from line center was readily evidenced by amplitude and 
shape changes in the laser pulse. The output of the laser passes 
through an amplifier tube (G) which provides a pulse of typically 
2kW peak power and 250 nsec width at half maximum, a dimethyl 
ether absorption cell (H) utilized to vary the pulse intensity, 
and two C0 2 sample cells (I,J). It is then attenuated by filters 
(K) and focused by a lens (L) onto a Ge:Au detector (M) used 
to monitor the laser intensity and pulse shape. The 4.3y 


fluorescence from the 00 °1 to 00 °0 CO^ level (See Fig. 2c) 
was observed with InSb detectors (N) . The second sample cell 
(J) , containing CO 2 at relatively high pressure (4 torr) , was 
used to insure that nonlinearities observed in the low pressure 
sample were not due to undetected change in the laser frequency 
or pulse shape. 

In order to obtain uniform intensity in the low pressure 
sample cell (I) , a 5mm dia. aperture was placed inside the cell 
in the center of the beam. This, together with a 7mm aperture 
placed 4 cm behind it, limits the field of view of the detector 
to this central 4 cm region, a distance in which diffraction of 
the laser beam does not produce serious nonuniformities. These 
apertures also served to reduce a linear background signal be- 
lieved to orginate from adsorbed on the cell windows. The 
laser profile at the position of the first aperture was measured 
to be uniform to within 30%, which gives a sufficiently small 
field variation to avoid severe smoothing of the experimental 
curve. Attenuation of the laser resulted in no detectable change 
in the laser profile. 
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The possibility of a frequency chirp during the laser puls 
which would greatly alter the nature of the observed effect (see 
Fig. lb) , was checked both by beating the laser against a cw 
CC >2 laser and by splitting the pulse and recombining it in an 
interferometric arrangement with a 12m delay in one arm. The 
latter technique is sensitive to small frequency changes during 
the pulse , provided the pulse shape is stable. By careful ad- 
justment of the position of the aperture near the rotating mirror 
the frequency chirp could be kept under 2MHz during the pulse 


(dw/dt<4MHz/visec) . 

Although fluorescent signals were directly observable on 


This sensitivity compares well with magnetic resonance. 

Experimental data (Fig. 3) reproducibly shows modulation of 
the fluorescence intensity at a pressure of 0.05 torr. This 
determines the laser intensity in terms of 0 . In order to compare 
our data with theory, the above analysis must be modified. In 
C0 2 , rotational relaxation is the main T^ mechanisms once a 
resonant molecule suffers a collision which changes its rotational 
state, it no longer interacts significantly with the radiation 
field. It nevertheless contributes to the measured signal since 
fluorescence is observed from all rotational levels of the 00°1 
vibrational state to the ground state. The proper expression 
for comparison with our data is thus 


One effect of this additional term is to cause the modulation to 
disappear more rapidly as the relaxation times are made shorter 
than Tp. Also, the fluorescence increases for shorter T 1 be- 
cause rotational thermalization tends to restore the population 
of both interacting levels to their equilibrium values. 

An additional complication is the presence of radiation 
trapping, which introduces a pressure dependence in the rela- 
tionship between observed fluorescence and excited population. 


an oscilloscope even at our lowest pressure (0.05 torr), signal 

14 

averaging was necessary in order to obtain accuracy. It is 
interesting to note that at 0.05 torr, in the region of y£Xp/ft = it 

9 

the signal originates from approximately 2 x 10 excited molecules 



( 6 ) 


The data plotted in Fig. 3 are proportional to the intensity 
of fluorescence immediately after the pulse, normalized to CC> 2 
pressure and multiplied at each pressure by a factor R(p) 
introduced to account for radiation trapping. At low pressures, 
a correction has also been made for a residual linear background 
signal. Theoretical curves are shown utilizing the known relaxation 
rate of 1/T 2 = 2 x 10 /torr-sec obtained from pressure broadening 
studies , ^ ' ^ 6 together with a value of 1/T^ = 10^/torr-sec^ in 
Eqs. (2) and (6). The set of numbers R(p) which produce 
agreement with these theoretical curves (see Fig. 3) is consistent 
with the trapped lifetimes measured by Kovacs and Javan. 6 
However, the uncertainty in this trapping factor for our cell 
geometry and the insensitivity to T^ of the curves resulting from 
Eq. (6) prevent an accurate determination of T^ . Variation of 
T 2 by 50% produces significantly worse agreement with the set of 

theoretical curves. The measured x p of 0.56 ysec and energy/unit 

2 ... 
area of 0.024 mj/cm in the region of the 2ir condition give a 

• -20 

value of y (P20) = (2.0 ±0.4) x 10 esu, in agreement with 

U /j 

the value of 1.7 x 10 -20 esu deduced from absorption measurements. 

The example discussed above shows the potential of the 
fluorescence method for the measurement of the dipole matrix 
elements and the relaxation parameters T^ and T 2 in gases in 
the presence of inhomogeneous broadening. Photon echoes, self- 
induced transparency and optical nutation measure similar parameters 
and thus should be compared with our method. When applicable the 
fluorescence method has the advantage of simplicity and sensitivity: 
it can measure relaxation in a sample of total absorption 10 • ** 
and relaxation times as short as the pulse length itself. This 
feature makes the method attractive for the measurement of short 


relaxation times, especially using tunable short-pulse lasers. 

Despite the advantages claimed above for this technique, 
we think it necessary to point out certain experimental difficulties. 
It is necessary to have a stable, reproducible, single mode laser 
with negligible frequency chirp and uniform amplitude across the 
beam. It is also important to eliminate any linear background 
signal arising, for example, from gas adsorbed on cell windows. 

On the positive side, we note that the magnitude of the dips should 
be much larger than those reported here if the condition x <Aw D ^ is 
fulfilled. 

Two simplifying ' assumptions were made in the derivation of 

Eq. 2: the third level has been neglected, and the radiative 

reaction on the level pair 1-2 has not been included. It can be shown, 

using the methods of Feld and Javan, ^ that under the influence of the 

strong field the radiative lifetime of the level pair 2-3 can change 

1 8 

only up to a factor of 2. Superradiant phenomena can also play an 

important role in determining the shape of the fluorescence curve as 

well as the fluorescence decay envelope, due to the enhanced rate 

18 1 

arising from the' superradiant state. f Both effects are negligible 
in our experiment. 

In conclusion, we note that this technique is also applicable to 
the study of relaxation in solids and in electronically excited 

19 

molecular states which are coupled to a quasi continuum of levels. 

It can be extended to include the time development of the angular 
distribution and polarization of the emitted radiation, which directly 
measures the equilibration of the- degenerate M-sublevels. 

We wish to thank Prof. A. Javan for his continuous encouragement 
and help in this work. The initial phases of this work were done in 
cooperation with Dr. C. K. Rhodes, who also bequeathed much of the 
apparatus . 
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FIGURE CAPTIONS 

I 

I 

Fig. 1 . Population of upper level as function of 0, calculated 

for various pulse shapes (S=sc]uare, H=approximate hyperbolic- 

secant, A=approximation of actual pulse of Fig. 2b) with 

t =0.6 ysec and for Aw_/2tt=60 MHz. (a) No relaxation. P,Q 
P ^ 

indicate average population of degenerate sublevels of high 
J, AJ=±1 or AJ=0 transition. (b) Effect of frequency chirp 
of 4 MHz/ysec. (c). Effect of relaxation, labeled by T 2 ,T X 
values (in ysec) as follows: B=°°,D=1, 1=0.1, L=0.0125. 

Fig. 2 . (a) Schematic of experimental apparatus (see text for 

details). (b) Multiple exposure of laser pulses, 100 nsec/div. 
(c) Relevant transitions in C0 2 . 

Fig. 3 . Observed fluorescence normalized to C0 2 pressure 
and multiplied by trapping factor R(p) as function of square 
root of laser intensity (in units of 0). Theoretical 
curves from Eq.6, labeled by T^^ values (in ysec) as 
follows: C=2, D=l, E=0.5, F=0.25, G=0.2, 1=0.1, J=0.05, 

K=0 . 025 , L=0 .0125. Origin is displaced for clarity. Error bars 
indicate typical uncertainty. 
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LASER- INDUCED LINE-NARROWING EFFECTS 
I N COUPLED DOPPLER-BROADENED TRANSITIONS . 
II: STANDING WAVE FEATURES 


B. J. Feldman and M. S. Feld 


Abstract 


Previous theoretical results on the influence of a laser on 
the line shape of a coupled transition — laser induced line nar- 
rowing -- have been restricted to the case where the laser is de- 
tuned from the center of its atomic gain profile or is in the form 
of a traveling wave. This paper extends those results to include 
the case where the laser is an intense standing wave field tunable 
to the center of its atomic gain profile (conditions for Lamb dip) . 
The effect may be observed either in transmission, by probing the 
coupled transition with a weak traveling-wave field coaxial with the 
laser field, or else in spontaneous emission from the coupled trans- 
ition viewed along the axis of the laser field. It is now well- 
for 

known that^a laser detuned from the center of its atomic gain pro- 
file, two narrov Lorentzian resonances of different widths appear 
superimposed upon the broad background signal at frequencies sym- 
metrically located about the corresponding line center. When the 
laser is tuned to the center of its gain profile, however, additional 
fine structure develops. This structure, which is particularly sig— 
n ificant when the laser field is intense, may have important appli- 
cations in high resolution spectroscopy and laser frequency stabili- 
zation. In this paper the laser frequency may be smaller or larger 
than the frequency of the coupled transition. In the latter case an 
intense laser introduces additional splitting effects, even when 
the laser is detuned. Splitting effects due to weakly saturating 
laser fields are also discussed. The problem is formulated by ex- 


panding elements of the ensemble-averaged density matrix in infinite 
series of spatial Fourier components. A perturbation technique is 
employed, valid for a weak probe field and a standing wave field 
of arbitrary intensity. One obtains an expression for emission induced 
by the probe field due to atoms moving with one velocity, written 
in terms of continued fractions in the general case and with 
Bessel functions in an important special case. This expression is 
integrated over the atomic velocity distribution by means of com- 
puter to obtain the total emission due to atoms moving with all 
velocities. A detailed discussion of line shape and of the physi- 
cal processes involved is included. 
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I. Introduction 

The lineshape of a Doppler-broadened transition is dramatically 
altered by the presence of a standing-wave laser field resonating 
with a second Doppler-broadened transition sharing a common level 
(Fig. la). For a laser detuned from the center of its atomic gain 
profile, two narrow Lorentzian resonances of different widths 
appear superimposed upon the broad background signal of the coupled 
transition at frequencies symmetrically located about its line 
center (Fig. lb).'*' 6 Recent observations of this effect, ^ 11 
called "laser induced line narrowing", confirm the predictions of 
the theory. The line-narrowing produced, which can be 100-1000 
times narrower than the Doppler background, has been utilized in 
a variety of ways as a high-resolution spectroscopic technique 
for determining isotope shifts, 12 fine 13 and hyperfine 11, 14 ' 15 
structure, g-factors 14, ^ and linewidth parameters. 1 ^' 1 ® / ^ 

The effect may be observed either in transmission, by probing the 
coupled transition with a weak traveling-wave field co-axial with 
the laser field, or else in spontaneous emission from the coupled 
transition viewed along the axis of the laser field (Fig. 2) . 

Note that identical lineshapes will be observed in both cases 
since the spontaneous emission lineshape is equivalent to the line- 
shape of emission induced by a probe field containing a single 

4 

photon in each mode. 

The overall features of the line-narrowing effect may be 
understood in terms of the velocity selection of moving atoms 
by the standing-wave laser field. The field selectively interacts 
with atoms whose velocities Doppler shift one of its traveling 
wave components into resonance. This produces changes in the 
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laser level populations over two narrow intervals symmetrically 
located about the center of the velocity distribution. These 
changes reflect themselves in the gain profile of the coupled 
transition. The above explanation does not account for the different 
widths of the two change signals, nor, as shall be seen, for the 
other striking lineshape features to be described below. The 
formulation, analysis and interpretation of these effects form 
the content of this paper. 

Previous theoretical results 1 6 on the influence of the laser 
on the lineshape of a coupled transition have been limited to the 
case where the laser is detuned from the center of its atomic gain 
profile (or is weak). A detailed treatment was given in Ref. 4, 
the first paper in this series (which shall be referred to here- 
after as "I"), using a method of calculation in which the E.M. 
fields are treated classically. The present paper extends this 
approach to the important case in which the laser is of arbitrary 
intensity and may be tuned to the center of its gain profile 

(conditions for Lamb dip) . This case is of particular interest 

20 21 

in view of the fact that recent theoretical studies ' predict 

that under these conditions the laser-induced change in the velocity 

distribution departs from a simple Lorentzian-like curve. Instead, 

. 22 

additional fine structure develops for atoms with low velocities 
(Fig. 3a) . This fine structure does not manifest itself in the 
Lamb dip. We shall show below that related fine structure does 
appear in the gain profile of the coupled transition. This 
structure is particularly significant when the laser field is 
intense. An example is shown in Fig. 3b, where the fine structure 
is clearly discernable on the wings of the central dip. The 


4 


additional central feature of Fig. 3a results from other spatial 
interference effects to be described below. 

When the laser field is intense other related lineshape 

effects can occur, even when the laser is detuned from the atomic 

center frequency. For example, when the laser frequency is higher 

than the frequency of the coupled transition under study, the 

narrow resonance splits into two (Fig. 4). This new effect should 

not be confused with another splitting effect which occurs for 
, 23-25 

weak fields. Both of these effects are treated below. 

It should also be noted that the lineshape details analyzed 
here are of importance in considering stabilization schemes based 
on laser induced line narrowing effects. 

This paper is divided into five sections. In Section II we 
set up the problem using the ensemble-averaged density matrix 
equations of motion. We obtain an expression for the small- 
signal gain profile of a transition coupled to a laser transition, 
valid for arbitrary values of laser intensity, laser detuning, 
level lifetimes, Doppler width, and relative frequencies of laser 
and coupled transitions. These expressions are given in terms 
of continued fractions in the general case and Bessel functions 
in important special cases. In Section III computer methods of 
evaluating these expressions are discussed. In Section IV we 
present the results in a number of special cases which emphasize 
the different physical effects which come into play. These results 
are compared with a simplified expression obtainable from a rate- 
equation approximation, and also with the independent field 
approximation . A brief discussion of the area (integrated intensity) 
properties of these results is given in Section V. Lengthy 


5 


mathematical discussions which would interrupt the continuity are 
deferred to appendices. 
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II. Theory 

The theoretical problem under consideration is the interaction of a 

coupled Doppler -broadened three-level system with two co-axial applied 

fields, an intense standing-wave field resonating with the "laser" transition 

and a weak traveling-wave probe field resonating with the coupled transition. 

To be concrete we shall deal with the folded level configuration 

of Fig. la. However, the results may be easily extended to the 

26 

cascade configuration. 

In formulating the problem it is most convenient to treat 
the applied fields classically. It is important to note, however, 

that the resulting line shape expression is also applicable to the 

£ . . 27 
case of spontaneous emission at the probe frequency. The formalism 

of the ensemble-averaged density matrix is adopted. In this 

formalism one considers the interaction of the applied fields 

with atoms moving with a given axial velocity. Having obtained 

the response at the probe frequency, it is necessary to sum up 

over the entire distribution of atomic velocities. Note that 

in the ensemble-averaged formalism the initial conditions at 

which individual atoms are produced have already been averaged 

over, a considerable simplification. 

The problem is solved by means of a perturbation technique. In the 
absence of the probe field the problem reduces to the interaction of a standing- 
wave field with a two-level system. We then consider the influence of the 
probe field as a small perturbation on the unperturbed system. The un- 
perturbed problem has already been solved for arbitrary intensities and 
detuning of the standing-wave field in Ref. 21. Thus, we may directly 
substitute the required terms into the perturbation equations. In the unperturbed 
solution it is important to include spatial variations of high harmonic content 
produced by the standing-wave field in the level populations and the 
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induced polarization. In Ref. 21 the unperturbed solution is 
expanded in an infinite series of spatial Fourier components. In 
the present problem it is also convenient to expand the perturbation 
solution in spatial Fourier series. This immediately leads to 
an infinite set of coupled nonhomogeneous linear difference equations 
which may be solved subject to the appropriate boundary conditions. 
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In the general case the solution may be expressed in terms of continued 
fractions and in important special cases in terms of Bessel functions. One 
then obtains an expression for the emission induced by the probe field due 
to atoms moving with one velocity. This expression is integrated over the 
atomic velocity distribution by means of computer to obtain the total 
emission due to atoms moving with all velocities. In certain special cases 
the integrations may be done analytically. 

The three level system under study is of the type shown in Fig. la. Level 

0, the common level, is coupled to levels 1 and 2 by electric dipole matrix 

elements (jl^ q and p 2Q , respectively. From parity considerations, p 2 i =0 and 

there are no diagonal matrix elements. Denote the energy of level j by 

fiW , and let W,..- W =w :u>, and o 0 fall in tne optical-infrared region. We 
3 0 3 3 1 ^ 

shall assume that jco 2 | is large compared to the Doppler widths so that 
E 1 does not resonate with gj 2 , nor E 2 with c^. The system interacts with a 
strong standing wave laser field E.^ (z,t) at , a frequency close to w 2 - 

The resonance at is probed by the weak field E^(z,t) at variable frequency 
. Specifically: 


E 1 (z,t) = fi° 1 cos ( Q jt-kjZ +<{> x ) , 


(la) 


o 


L 2 (z, t)-E 2 sink 2 z cos (^ 2 t "t"(jx 2 ). 


(lb) 


with 



(lc) 


and <p x and <p^a.re constant phase factors. The additional condition 



(Id) 
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(L-cavity length, m is an integer) is imposed on E 2 by the laser resonator. 
The total Hamiltonian for the system is 


H=H+V. 


(2) 


where H Q is the Hamiltonian of an unperturbed atom with stationary states of 


energy f ^ , and V is the interaction Hamiltonian, 

j 

V = -jiE(z.t), 

E(z,t) = E^(z,t) + E 2 (z,t), 
with electric dipole operator p. 

The problem is most readily treated by means of the ensemble 
averaged density matrix p(v, z,t), which describes an ensemble of atoms 
at coordinates {z,t} moving with axial (z-axis) velocity v. The equation 
of motion of p, derived in Ref. 21, is given by 


(3a) 

(3b) 


(rr +v P - !/2 {p - p°, r\ 


(4) 


where [] and {} are commutator and anti-commutator brackets, respectively 
The left-hand side is the total derrivative of o . Relaxation is accounted for 
in the second term on the rijit-iand side through r, which has elements 

Wv (5) 

O 

where yds the relaxation rate of level j. Also,P describes the steady- 
state ensemble, assumed spatially uniform, in the absence of applied fields: 


<V n i (v) V 


(M 


n.(v) being the steady-state population density per velocity interval in level 

J 

j in the absence of applied fields. We may write 


n j'v> - N J G / V) ' 

with N. the total population density of level j and G (v) its normalized 
J J 

velocity distribution: 

«• 

■( v)dv = 1. 

J 


(7a) 


(7b) 
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For further details regarding Eq. (4) see Ref. 21. 

The partial induced polarization due to the atoms moving with 
velocities in the interval between v and v + dv is 

P(v, z,t)dv = Tr [>iP]dv. (8) 

The net polarization is obtained by considering contributions from all 
velocities: 

00 

P(z,t) Tr [ji 0 ]dv. (9) 

— 00 

To formulate the perturbation approach let us rewrite the 
Hamiltonian of Eq. (2) in the form 

H = H U +h' , (10) 

with 

H U = H -|iE 9 (11) 

o z 

and 

H' = -pE r (12) 

Thus, H u is the unperturbed Hamiltonian, describing the interaction of 
the standing wave field with the atomic system in the absence of E^. The 
corresponding unperturbed density matrix, p u , is the solution to Eq. (4) 

I 

with H equal to zero: 

(■gr +v a4 > P u= - T p u ] - !/ 2 r}. (13) 

f 

The additional term H , due to the weak probe field, is taken to be a 
small perturbation on H . Thus, H produces a small correction p to p u ; 
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p 



(14) 


Inserting Eqs. (10) and (14) into Eq. (4), subtracting Eq. (13), and 

i 

neglecting the [H , p'] term, which does not contribute to p'to lowest 

I 

order in H , one obtains 

t 3t +v > »' -*• l H ' >° U 1 - im u . p' ] - >/ 2 fe* »r>. (15) 

Assuming p U is known, p' is completely specified by Eq. (15). 

Equation (13) has been studied in detail in Ref. 21. Its solution, p u 
valid for arbitrary values of E. ? , has elements of the form 


u 

°20 


(v, z,t) 


ii 1 ft, at 

= X U (v, z)e 


(16a) 


X U (v, z) 


n 


02 


7 [ v 

o3a“n>l 


ink_z -ink-z 

e ^ +lf e 2 
n 


]• 


(16b) 


u u 
p 21 = p 01 = 


0, 


(16c) 


u _ u* 
p ik ~ p ki’ 


(16d) 


Zb 

even n>0 


ink£Z 

(v)e +c. c 




p 00 (v ' M )= n 0 + n 0 2 


(I6e) 
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P22 (v ' M) = n 2 + n 02^[* b n (V)eink2Z+C ' C ‘ J ' 

even n>0 


u 

PH n i) 


where 


n. . = n. -n . • 
ij i J 


(16f) 


(I6g) 


(16h) 


(In obtaining Eqs. (16) it is assumed, as stated above, that E 2 does not 

resonate withw .) For convenience a , b andH^- , obtained in Ref. 21, 

1 ' n’ n n 

are given in Table 1. Note that Z n Eq. (T 14) , is a continued 
fraction which is to be evaluated subject to the boundary condition 
Z n ->-0 for n-*- 00 . The computation of Z^ is discussed in Section III. 

Let us now consider Eq. (15). The weak probe field does not perturb 
the level populations, so that 

Pkk = °* k= °’ 1 ’ 2 - (17) 


Therefore, using Eqs. (16c) and (I6g), the equations of motion of the off- 
diagonal martix elements of p may be written as 

("it + V 8~z + la) l +7 lo) p 01 = "T J^01 E l(- D 00" n l) _M, 02 E 2 P 21 j’ (18a) 

(^"t +V %. +l(a) r u 2^ +7 2l) p 21 = " ~t [^01 S 1 P 20 " M '20 E 2 p 0lJ > (18b) 

f^t +V 5"z ‘ lt0 2 + Y 2oj p 20 = ~~* ^20 E 2^22' p Oo) + ^ 10 E !p 2 i]- (18c) 


This set of equations may be solved by inserting E^ and E 2 [Eq. (1) ] in 
complex form: 
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EjU.t) = Aje-Ve^l'tc.o, 


(19a) 


E 2 (z,t) = (A 2 /2)(e lk 2 z -e‘ 1k 2 z )e 1 ' n 2 t + c.c. 


(19b) 


with 


A -ApO 

A 1 " iE i e , 


(20a) 


*2-T& 1 ** 


(20b) 


The various spatial and temporal Fourier components of E^ and E 1? drive 

the off-diagonal elements of p Let us consider the variations of the o' -j. 

Of special importance are the coefficients of p' . - on the left hand sides of 

J 

Eqs. (18). These are associated with the resonant behavior of the induced 

polarization. The only important frequency components of a o' . • are those 

3 

for which 9p'^j/3t cancels the iw factor for appropriate values and 
, thereby reducing the coefficient of p'_ toy.y-va/az . Therefore, a 
nearly complete solution to Eqs. (18) will be of the form 


°oi = A(v - z )e_1 ^ 


(21a) 


' = D(v, z)e i(fi 2 _fi l )t . 


P 21 


(21b) 


’20 


\(v, z)e’ ^2*, 


(21c) 


where D and X. are time -independent. Inserting Eqs. (16a), (19) and 
(21) into Eqs. (18) and equating like coefficients of the 

complex frequency factors, one finds that*- (v, z) = 0 and that 
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(A l +iv ii +i no )A=a * eiklZ( Poo' n l ) + ji3(e ik 2 Z -e" ik 2 z )D, (22a) 

(A 12 +iv 3 - + i 7 21 )D= Qf*e lk l Z X u +-1 3*(e lk 2 Z -e _lk 2 z )A, (22b) 

where 


a " ^10 A ]/ ^ , 

(23) 

* 


3 = " 1 ^02 A 2 // ^ • 

(24) 

A .= ft .-co • 

3 3 3, 

(25) 

A 12 = A 1 ‘ A 2* 

(26) 


We may proceed by expressing A and D in Fourier expansions in k 2 z. Note 
in Eq. (16) that Pqq is expanded in the even harmonics of k 2 z, whereas \ 
is expanded in the odd harmonics. It directly follows from Eq. (22) that 
A will consist of the even harmonics of k 2 z, and D of the odd harmonics. 
Therefore, the required solution will be of the form 


a / \ ^ i k, z 

A( v, z)=a e 1 


A 0+ £{A>>e‘"V. A >, e ""V 1 


even 

n>2 


D(v, z)=a 




(v)e ink 2 2 +Di ;(v)e- ink 2 z 


odd 

n>l 


}• 


(27a) 


(27b) 


The induced polarization at ft ^obtained from Eqs. (8) and (21), is 
of the form 


p^v, z,t) = Re^2p 10 


Ae' ifi l 


*] 


(28) 
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The average power increase per unit volume at induced by Ep due to 
atoms located at z and moving with axial velocity v is given by 

AI(v,z) dvdz = -^(v.z.tJE^z.tX dvdz. 

time average 

The net power increase per unit volume is obtained by integrating this 
quantity over the distribution of atomic velocities and the length, l of the 
active medium/ centered at position a: 


a+1/2 


A I = 


dz 


a.- 1/2 J'oo 


J' dvAI(v, z). (30 a) 

— — , _ CO 

Making use of Eqs. (29) and (27a), and indicating velocity averages by / • i.e. 

00 

< 0 > = f 0 (v) dv, 

i. < A± n>l ' 

(30b) 


— oo 

we obtain 

A I =-2fi& 


l l a I ^ 


Im <A _> + additional terms of the form 

0 nfc 


where X=2ir/k^ is the wavelength of the intense standing-wave field. Since 

. + 

\ / 1 £< 1 , the additional terms are negligible as long as the do 

not diverge for large n. As shown in Appendix A , the cave rge nee of the 

A* as n— 9 oo is a necessary condition for a physically acceptable solution. 

Accordingly, the additional terms may be neglected, and we obtain 


AI =-2 fiQ 1 1 ct| 2 J£Tm<A 0 >. (32) 

Note that only Aq contributes to the net emitted power. 

We now solve for Aq . Inserting Eqs. (16b), (I6e) and (27) into Eqs. 
(22) and equating like harmonics of k^ z , we obtain a set of equations which 
can be written in the following compact form: 


£x + (n+l) - x + (n-l) J + L(+n)x ± 


(n) = n 02 t+(n), n>T, 


(33a) 


i-^JV^l) - x_(l) ] 


X (1) - x_(l) + L(0)x{0) = t (0). 


(33b) 
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The following definitions have been made: 


f^//3o> D n 


x + (n) =< 


X + (0) = x(0) = A q 

n + n 


t + (n) = 


V 


n 




t (n) =< 


a* 

n 


n = 1, 3, 5. 
n = 2, 4, 6. 


n = 1, 3, 5. 
n = 2,4,6. 


r. = 1. 3. 5. 


n = 2, 4, 6. 


(34a) 

(34b) 

(34c) 

(35a) 

(35b) 


(35c) 

(35d) 


t(0) n Q2 [a Q + a£] + n Q1 


f 

L(+ n) =< 


^12 ” ^ k l ink 2^ v + ^21 n = i* 3 * 5 --- 
- (I^+ni^) v + i Tl0 n = 0, 2,4,6. 


(35e) 


(36a) 


(36b) 


Equations (33a) and (33b) constitute an infinite set of coupled nonhomogeneous 
linear difference equations, in which the t + (n)' s , known quantities 
determined by the unperturbed solution, act as driving terms for the x + (n)' s. 
We are particularly interested in obtaining A Q since, from Eq. (32), 
the emitted power is proportional to ImAg. It is shown in Appendix A 
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that the solution for Aq is 


n 0l 4 ' n 02 


o - -JU 


^o +a *o + 

r u + (l)-u.(l) J 


+ L (0) 


( 37 ) 


where 


U ± ( J) = K 1 w + (k)> 


(38) 


and W (k) is the continued fraction 
1 ±^0 1 


W + (k) = 


2 L i. k] 


A 


1- 


L[±k] 




a 


1- 4 L[+(k+I)1 L'F(^+2y 


1 - 


(39a) 


1 - 


evaluated subject to the condition 


lim W± (k) = 0 . 

k-» oo 


For the special case 7 0 = T 2 / A -=O we find 


(39b) 


u ± (j) . , 


J t (0) ( - S o /k 2 v) 


(40a) 


t ( j ) = +L (+j )/k 2 v, 


(40b) 
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where J ^z) represents a Bessel function of complex order t and (real) 
argument z. The final expression for AI may be obtained from Eqs. (32), 
(37) and (38). The numerical evaluation of W+(k)i S discussed in Section III. 
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III. Computer Evaluation of AI 

The calculation of AI , the emitted power at ft p from Eqs. (35), 
(37), (38), (39) and ( 31)/requires the evaluation of continued fractions and 
velocity integrals. The computations have been performed on the M. I. T. 

360 computer using double precision complex arithmetic. Continued 
fractions were evaluated by comparing the value obtained by truncating the 
continued fraction after j terms with the value obtained by truncating after 
j+1 terms. The algorithm used in calculating a truncated continued fraction 
of j terms is given in Appendix C. If the absolute value of the fractional 
differen ce between the two calculated values is less than 10 the j+1^* 1 
value is taken as the numerical answer. The u + 's and y' s represented as 
convergent continued fractions have been calculated from Eqs. (38) and Table 
I with all the continued fractions computed independently to insure that small 
errors did not propagate. In the summations in Eqs. (37) sufficient 
terms were included so that the absolute value of the ratio of the last value 
to the first is less than 10 The convergence of this sum is a consequence 
of the boundary conditions, which require to approach 0 for large n 
(Appendix A). 

As a check, when anc * ^2 = ^1 we have compared our 

continued fraction results with the Bessel function expressions and find 
agreement to 6 places of accuracy. 

To simplify the interpretation of our results we have chosen a 
rectangular population distribution. 


for all values of j. This procedure greatly reduces the effects of a finite 
velocity distribution, which would arise had we chosen a more realistic 




( 40 ) 
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2 8 

population distribution. It would, of course, be easy to insert 
any desired form of (v) into the program. 

The velocity integrations are evaluated using Simpson's 
rule with intervals of .2k2v/Y2Q. The accuracy of the integration 
is better than 2%. The accuracy of our calculations was checked 
by decreasing the step size and increasing the limits of integration. 
Our numerical results are in excellent agreement with the lim- 
iting case of a weakly saturating standing wave field (3rd order 
polarization) and with the case in which the laser is of arbitrary 
intensity but significantly detuned, cases for which we have 
obtained analytic expressions. (See discussion below.) 
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IV. Results 

Equation (32) may be written in the following form: 
AI ( %) = &ft i<( n 2 - n i> J 21 + (n 0 - n i> J oi> ' 

(as before, <> indicates velocity average), where 

J 21 =2ja | 2 ilm 1 • < L 


£1 


J 01 =-2 | a | 2 £Im 


i^jv i)_u - (1) ] + L(0) 

f 1+a o +a o + + i<3) u .(j) J 


^ £u + (l)-u_(l) j + L(0) 


(41) 


(42a) 


(42b) 


Expressions in the same form for the off-resonance case 
1/2 

(A 2 /y 2 q ) > d+l/2) were previously discussed in I, where it 

was shown that represents the double - quantum transition rate 

per atom, in which atoms produced in level 2 subsequently decay 

from level 1 by the net exchange of two photons with the applied 

radiation fields, one absorbed at ft_ and one emitted at ft, .J was 

2 1 01 

introduced as the single - quantum transition rate per atom, due 
to atoms produced in level 0 which subsequently decay from level 1, 
emitting a single photon at fi^. In this paper , however , we find 
it more convenient to express our results for Al(ft^) in terms of 
ng-n 2 and n^-n^. We thus have, equivalent to Eq. (41), 
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) * A 1 £^ R 02? + ^ R 01^~| ' 


where 


R 02 = ‘ 2n oi arita 


laf« 


vv 


oo 

^(1) J j~t + ( j) u + (j) +t_(j)u_(j) J 
^u + (l)-u.(l) j+L(0) 


R oi = - 2n oi°< 1 2 * * ^ <— T 


[u + (l)-u.(l)J + L(0) 


(43a) 


(43b) 


(43c) 


In Figs. 5 through 8, <R o 2 > and <R 01 > are P lotted for various 

parameters. For convenience we have introduced the usual saturation 

parameter 2 

2B o 

1 = Y Y ' (44) 

0 2 

as a measure of the intensity of the standing wave field 
(E^). In Figs. 5/ 7a and 7b the standing wave field is 
on resonance In F ^9 S * 6,7c and 7d the standing wave 

field is off resonance 02 >> ( l +I /2 ) 1//2 ] • In Figs. 8a 

1/2 

and 8b the standing wave field is near resonance [A 2 /Yq 2 ~ (1 +i /2 ) ]• 
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Brief ly, the general behavior of the curves is as follows: 


one broad, one narrow, but of equal area, occurs symmetrically 
located about In the folded configuration discussed here 

(Fig. la) the narrow resonance and the laser are always detuned- 
to the same side of their respective Doppler profiles. (In a 
cascade configuration the narrow and broad resonances would be 
interchanged [Ref. I]). In the special case k 2 >k^ narrow 

resonance splits into two for intense laser fields. As the laser 
is tuned to line center the resonances overlaps and additional 
complex structure develops. 


laser when detuned, except in the case in which structure 

develops on the same side of the Doppler profile as the narrow 


on the opposite side of the Doppler profile.) As the laser is 
tuned to line center structure emerges for strong laser fields. 
In all cases, the area under the background curve is unchanged 
by the presence of the laser field. 

The additional features occurring on resonance are due to 
the standing wave nature of the intense field, whose traveling 
wave components simultaneously interact with atoms moving in the 
same narrow range of velocities. It will be shown below that 
off resonance, where the traveling wave components of no 
longer couple to the same atomic velocity bandjthat the line 
shape can be analyzed in terms of the two oppositely directed 
traveling waves at ^ which do not couple to one another. 



curves : 


When the laser is detuned a pair of resonances, 



(background) : The background is unaffected by the 


resonance in the ( r q 2^ curves . (In no case does structure occur 
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Intense Standing Wave Field - On Resonance [A ^ -0 ] 

Typical results obtained from Eq. (43) are shown in Figs. 

5a-5f. To understand the significance of these results and 
how they differ from the curves produced by an intense traveling 
wave field, consider the case k^=]< 2 . For the traveling wave 
case broad (E^ j | ) or narrow (E^ anti-| j E^ ) resonances are 
produced in the curves , depending on the relative directions 
of the two fields. No additional features occur as the traveling 
wave laser field is tuned through resonance. In the standing 
wave case the broad and narrow resonances occur simultaneously 
on opposite sides of the gain profile when the laser is 
detuned. When the laser approaches line center, these resonances 
overlap. The resulting curve is not simply a superposition of 
broad and narrow signals. Complicated structure occurs. This 
structure results from the fact that the induced polarizations 
at associated with each of these resonances, being coherent 
with one another ^ interfere by means of spatial variations induced 
in both the population distribution [p^ Q Eq. (16e) ] and the 
parametric driving terms at ^2~^l tP r °duced by P 20 ' E< ^* ' 

causing complex line structure to develop [compare Eqs . (21a, b) 

and (27a, b) ] . 

On resonance the standing wave background curves also exhibit 
complex structure resulting from coupling of the spatial variations 
in the polarization at directly through p^, which drives the 
system parametrically at ^2 - ^l* Note that the characteristic 
behavior of the background term is distinct from that of the 
<^Rq 2 ^ term, occurring even when there is complete transparency 
at the laser transition (n 2 =nQ). It is also interesting to note 


-25- 


that the area under the background curves remains unaffected by 
the laser field in all cases, consequence of the fact that the 
area of the response of each atomic velocity ensemble is unaltered 
by the laser field. A comprehensive discussion of this fact 
will be presented in a forthcoming publication. 

It is worthwhile to reexamine the above line shape features 
from another point of view. In Refs. 20 and 21 it was shown 
that when an intense standing wave field is tuned to the center 
of the Doppler profile of a twc level system, fine structure occurs 
in the spatially averaged population distribution . An example 
is shown in Fig. 3a. This structure results from the coherent 
ringing of slow atoms moving through the spatial nodes of the 
standing wave field (See Ref. 21, Sec. 7 for detailed discussion.) 
This curve is compared with a typical <^R 2Q ^> curve in Fig. 3a. 

The influence of the velocity deformations in the spatially 
averaged population are clearly evident on the wings of the central 
dip in ^ r q 2 ^ • No such fine structure appears in the central 
tuning dip of the two level gain profile (or the Lamb Dip) , 
obtained by tuning the intense field through the center of the 
atomic resonance. Thus experimental observation of curves such 
as^Rg 2 ^>°f Fig. 3b would not only confirm the present results 
but would also provide evidence of the population structure 
predicted in Refs. 20 and 21. 

Indeed, a major incentive of this work has been to determine 
the extent to which the structure in the population distribution 
would manifest itself in the radiation at the coupled transition. 

This behavior was anticipated on the basis of a simple argument 
ignoring field interactions between and [Independent Field 
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Approximation (IFA) ] . In this limit the spontaneous emission 
at is given by 


Eq. (16e) , in the presence of the standing-wave laser field. 

One would expect the fine structure to wash out for Y^q >> Y 2 q 
since in this case the broad Lorentzian cannot follow the narrow 
structure within the population curve (whose fine structure is 
of the order of y 2 q). However, with Y^q£Y 2 q fi- ne structure would 
occur. Needless to say the IFA expression is not rigorous. For 
purposes of comparison the IFA results are plotted together with 
the exact expression for ^Rq 2 ^ [Eq • (43b)] and the corresponding 

rate equation approximation (REA), discussed below, in Fig. 9. 

For the case k^k 2 additional line shape features develop. 
These features are also observed when the laser is detuned from 
resonance, where they can be separated from the standing wave 
interactions discussed above. We defer discussion of these 
traveling-wave effects to the next sub-section, where algebraic 
expressions for ^Rq , valid for a detuned standing- 

wave laser field of arbitrary intensity, are presented. 

Intense Standing-Wave Field-Off Resonance (A 2 /v 2 q>> d+I/2) 1 ^ 2 ] 

Next consider the case when the laser field is off resonance 
(Figs. 6, 7c and 7d) . In this limit R Q2 and R Ql reduce to 



> 


( 45 ) 


where n Q +n Q 2 (a^ta^ is the spatially uniform component of p Q0 (v) 
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R 02 = ~ (1 / 2 ) n Q2^o I a I x 



R 01 “ 2n 01 


r 


|a| 2 *Im J ' 

l 



+ 



(46b) 


Where is given in Table I, Eq. (T21) . 

These results are equivalent to the Rate Equation Approximation 
(REA) , in which the spatial harmonics induced by the standing-wave field 
are neglected. (See Appendix A and Table I.) Since we are restricting 
our considerations here to the detuned case, Eqs. (46) further simplify 
to 


R 02 =R 02 (+ > + E 02 ( -> 


R 01 =R 01* R 01 <+) + R oi ( ->' 


(47a) 

(47b) 


with 
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R 02 ^ I 


‘02 



(48a) 


where 


1 _ 

F 


2 iv,„ + i 


Y 20 +i<4 2 lk 2 V) 


(48b) 


and 


R 01 (±)=-(l/2)n 01 |a| 2 6^Im 


[ 


L (0 ) I L(0)L(±1)- ^ 


L 


'01 01 


(49a) 


|ci| 2 llm 

r, I 



1 

[_ L (0 ) J 


(49b) 
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The terms naturally break into pairs (+ and -) as a consequence 
of the fact that when the standing-wave laser field is detuned 
from resonance the oppositely directed traveling-wave components 
of E 2 couple to distinct bands of atoms moving with opposite 

velocities. In fact, the (+)[(-)] refers to the component of E ? anti- 
parallel [parallel] to E^. Note that R^represents the usual background 

term in the absence of . These results were previously obtained 
in I, where the intense field was assumed to be in the form of 
a travelling wave. In the Doppler broadened limit the velocity 
integration of Egs. (48) and (49) may be carried out analytically. 


The results for the two cases k ± >k 2 and k^^ differ qualitatively 


Case 1: 


In this case we find that 


<R 01 >= R ° w oi<V k i K 


01 *1 


(50a) 




(50b) 


(50c) 


where 



(51a) 





(51b) 



(51c) 


(51d) 


= A 1 + (k 1 /k 2 ) A 2 


( 5?s) 


r 



(5lf) 


and 



(51g) 


Equations (50), valid for large values of E 2 , were first given 


with k 1 >_k 2 are plotted from Eqs . (43 b, c) in Figs. 6a,b,e,f,7c 

and 7d. It is found that the results agree with the 
analytical expressions, Eqs. (50). The <Rq 2 > curv es exhibit 
a pair of Lorentzian resonances of equal area at f2 + = 0 , of 
widths r + as in Eg. (50c) . Note that T + are significantly 
power broadened for intense laser fields ( Q >> 1 ) , as in Eq. (51f ) . 

The + and - resonances correspond to the interaction of E ± with that 
component of e 2 anti-parallel and parallel to E^, respectively. The 
<R 01 > curve exhibits no structure, in accord with Eq. (50b), 
and is identical to the background one would observe if E 2 


in 


I 


Curves of <r q^ > end <r q 9 > for the detuned case 


were absent. 
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Case 2: < 1^2 


In this case we find that 


^ R oi^ ~ ^ R oi + R oi^~^ “ R W 01 ( V k l> Im ( y ) ' 


(52a) 


< R 01 (+ )> = 0 , 


(52b) 


^ R 02 = Same as for k^ k 2 , Eq. (50c) , 


(52c) 


< R 02 ( - } > 


R Q n 02 (0) ( A /k ) B ° /*A 2 x 

R n Q1 (0) W 02 (A l /k l ) 2Q [k 2 j 


Im 


{ Y ' X (^ lr?~ 


i+2_ + < kl ' k 

2 T 


£r "* * *3 




-O k l~ k 2 

Q + i -~r- (A-+ 6) 
Y 0 K 1 


)} 


(52d) 


where 


l^i ■ 


(53a) 


T = ft_ + y_ + y 


1 , 


Y = + 


js 2 - T 2 j 


1/2 


(53b) 


(53c) 


-k. 


A + “ 2 (k 2 -k 1 ) (T - lY) + ( A 1 + iY 0l) 


(53d) 
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(53e) 


6 ■ k^q ( A 21- iY 2l) • 


(53f ) 


As an example, curves of 


<k 0 ,\, and (n oi y 


for the case k^=k^/2 are 


plotted from Eqs. (43b,c) in Figs. 6c, d and 7c, d. The results agree 
with curves plotted for the analytic expressions, Eq. (52), and 


curves still exhibit a Lorentzian resonance at ft + = 0 of width 
r + corresponding to the interaction of with the anti- 
parallel component of E ? , just as for the case k^> k 2 . However, 
the resonance at = 0, corresponding to the interaction of E^ 
with the parallel component of E ^ , exhibits a splitting which 
increases with increasing field strength E 2 and vanishes in the 
limit of weakly saturating fields, where the resonance reduces to 
a Lorentzian of width T_ (see following paragraph) . In addition 
to the features in structure in < R 01> also manifests itself. 

We see that a distinctly non-Lorentzian resonance appears at 
ft_= 0. This resonance features a central dip and broad wings 
peaking off center, in such a manner as to leave the area of the 
< R oi> background curve unaltered by the presence of E^ . The 
resonance is a result of the interaction of E^ with that component 
of E^ parallel to E^ , and as such is a travelling wave effect. 

This effect, which occurs even for weakly saturating values of E 2 , 
has recently been observed in Ref. 9. The splitting of the ^ r q 2 ^ 
resonance, which appears only for a strongly saturating value of 
E^ , is a new effect and has not yet been observed. 


differ significantly from those for the case k^>k£ 
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Intense Standing-Wave Field— Near Resonance L A 2^ Y 20 ~ (1 + I /2) 1/2 _ 

In the near resonance case interference effects produce 
additional lineshape alterations. Examples, plotted from 
Eqs . (43b, c) , are given in Figs. 8a, b. Note that these curves 

combine features of the off-resonance and on-resonance cases, 
leading to complex asymmetric lineshapes. 


-34- 


Limit of Weakly Saturating Laser Field (Q-l) 

The results for the case of weakly saturating laser are 
plotted in Fig. 7. The analytical expressions for the case 
k^>k 2 are the same in Eqs . (50) but with Q~ 1 ( indicating 

that the Lorentzian resonances of different width but of equal 
area occur in < ^Q 2 > sven in this limit. It follows from Eqs. 
(43b) and (43c) that in the weak saturation limit Eqs. (50) 
are valid even when the laser is tuned to the iine center. 

In this limit it follows from Eqs. (52) that the analytical 
results for the case k^<k 2 are identical to the results for 


^1-^2 t ^ ie exce P t; i- on <R 01 <-> > /Which is given by 


(23-25) 



, k 2 >k r (54) 


-35- 


Therefore, even for a weakly saturating laser field structure 
appears on the background lineshape when k^>k^ , although the 


appears when k^>k£ • Note also that the splitting of the narrow 


fields for ^2 > ^i ' disappears i n this limit. All these remarks 
are contingent upon the assumption that the Doppler widths are 
much greater than the natural widths . 

Computer plots of Eq. (43) for 1= 0.1 are given in 
Figs. 7c, d for the off-resonance case and Figs. 7a,b for the on- 
resonance case. Note that the latter curves are just the super- 
position at A 1 = 0 of the detuned resonances, indicating that 
spatial interference effects between components of the induced 
polarization at are not significant here. 


area under remains independent of E 2 . No such structure 



which occurs at intense laser 
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V. Equal Area Property 


Some interesting properties of the laser induced change 
signals emerge from the results presented here. First, as stated 
above, the area under the <R^> background is always independent 
of E^. Secondly, the area under the <R £ 2 > lineshape is 
independent of the field interactions between E^ and From 

this consideration it follows that when the laser is detuned 
the broad and narrow resonances are always equal in area, even 
when the narrow resonance splits in two. In fact, it may be 
shown that in each case the area under the exact lineshape curve 
{Eq. (43a)] is identically equal to that under the corresponding 
IFA curve [Eq. (45)]. This property is a manifestation of the 
fact, first pointed out by Javan^^ , that the area under the 
response of each individual atomic velocity ensemble is equal to 
that obtained from the corresponding IFA, even though the actual 
frequency behavior differs significantly. These results can be 
shown to follow from an elementary quantum mechanical analysis 
of the transition rates of three level systems coupled to applied 
fields. A detailed discussion will be deferred to a later 
publication. 
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Appendix A: Solution to the Nonhomogeneous Linear Difference 

Equations (33 ) 

Consider the coupled nonhomogeneous iinear difference equations, 
Eqs. (33) of the text, which may be written in the form 


x ± (n+l) -x ± (n-l) + A ± (n) x ± (n) = <J> ± (n), n >1 


(Ala) 


x + (l) -x_ (1) +A ( 0 ) x ( 0 ) = tf) ( 0 ) 


(Alb) 


x + (0) = x_ (0) = x (0) , 


(Ale) 


where 

A±(n) = 4^— L (±n) , ‘ (A2a) 


±2n 0 2 

4> ± (n) = — Ig — t ± (n) , 
p o 

A(0) = IT" L(0) 


and 

<t> (0 ) = j|-t (0) . 

6 o 


(A2b) 


(A2c) 


(A2d) 


Note that adjacent x + (n)'s are coupled together in three's. The 
<j> + (n)'s are prescribed driving terms. Equations (Alb) and (Ale) 
serve to join the x + (n) with the x_ (n) . Henceforth, the ± 
subscripts will be deleted when convenient. 


-39- 


The solution to Eqs. (Al) may be formed from the solutions to 

pp 

the associated homogeneous equations, x (n) . The homogeneous 
equations are obtained by setting the <J>(n) = 0 : 

x H (n + 1) -x H (n - 1) +A(n) x H (n) = 0, n >1 . (A3) 


Solutions to equations of the type Eq. (A3) are discussed in detail 
in Section 5 of Ref. 21. A brief account is given here in 
Appendix B, where it is shown that the most general solution for 

pp 

x (n) may be written in the form 


x H (n) = Cu(n) + Dv(n), 


(A4) 


where u(n) and v(n) separately satisfy Eq. (A3), and C and D are 
arbitrary constants independent of n. The u(n)'s and v(n)'s are distinct 
classes of solutions characterized by their behavior at large n. The 
u (n) ' s , which we call B-tvpe solutions, have the limiting behavior 


u (n) ~ K -(-& — Y 1 / n large. 
v ' nT^2ik 2 v J 


(A5a) 


The v (n) ' s , called N-type solutions, have the limiting behavior 

« 

v(n)~ (n - 1) I (- 2k ^~) * n large. (A5b) 


Accordingly, in the limit , u(n) vanishes whereas v(n) diverges. 

In subsequent discussions it will be convenient to normalize u(n) 
and v(n) such that 


u (0) = v (0) = 1. 


(A6) 
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Solutions to the nonhomogeneous equations, Eq. (Ala), may be 

31 

constructed using the method of variation of parameters. Let 
us assume x(n) to be of the form 

x(n) = C(n) u(n) + D (n) v(n), n >0 . (A7) 

Inserting Eq. (A7) into Eq. (Al), and eliminating terms containing 
A(n) by means of Eq. (A3), one obtains, for n^_l / 

6C(n)u(n + 1) + 6D(n)v(n + 1) 

+ 6C (n - 1) u (n - 1) + 5D(n - 1) v (n - 1) = 4> (n) , (A8) 

with 

6C(n) = C(n + 1) -C(n), (A9a) 

6D(n) = D (n + 1) -D(n) . (A9b) 

Trial solution (A7) does not uniquely specify the C(n)'s and 
D(n)'s, and in Eq. (A8) we are free to choose 

6C(n) u(n) + 6D(n) v(n) = 0, . n >0 , (A10) 

whereupon for n>JL, 


6C(n) u (n + 1) + 6D(n) v (n + 1) = <f»(n)* 


(All) 
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Combining Eqs. (A10) and (All), we find 


- 4> (n)v(n) 

° C(n) u (n+1) v (n)-u(n) v(n+l) 


n >1 . 


6D 


(n) = -<H n ) u < n > 

' n /n4-l ^ 17 / n ^ — 


u(n+l) v(n) - u(n) v(n+l) 


J 


(A12a) 


(Al2b) 


These expressions may be simplified by means of the following identity, 
obtained from Eq. (A3): 

u (n+1) v (n) -v (n+1) u (n) =- [u (n) v (n-1) -v (n) u (n-1) ] (A13) 

= (-l) n [u (1) - v(l)] . 


(The last step follows by induction. ) Thus, 


6C (n) 


(-l) n <t> (n) v (n) 
u(l) - v(l) 


6D(n) 


(— 1 ) n d> (n) u(n) 
u(l) - v (1) 



(Al4a) 


(Al4b) 


Accordingly, for n>l. 


u(l) - v(l) 


C (n) = C(l) + 


n-1 

E 

j=l 


(-l) 3 <Hj) V(j), 


(Al5a) 
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D (n) - DU) u(l)-v(l7 


n-1 

l 

j=l 


(-D^(j)u(j) • 


(Al5b) 


An additional equation, obtainable from Eq. (A10) for n = 0, is 


C(0) + D (0 ) = C(l) + D (1 ) 


(A16) 


since, by Eq. (A6), 


u (0) = v (0) = 1. 


(A17) 


Therefore, using Eq. (A6), we can write down the sequence 


x (0) = C(l)+D(l) , 


x(l) = C(l)u(l) +D (1 ) v (1), 


(Al8a) 


(A18b) 


{■ 


n-1 


x (n) = «|C(1) + u (1) -v ~ (1 T 4> ( j ) v (j ) (n) 


} 


+ j* D(1) " u (l) -vTTT ^ =1 


i n-l -i i 

1 E (-l) ] <Mj)u(j) Y v(n) ,r,n*l. (Al8c) 


} 


Thus the x(n) may be expressed in terms of two quantities independent of n, 

C(l) and D(l), as of yet undetermined. 

The value of D(l) may be established by considering the behavior of 
x(n) as . Recall from Section II that the x(n)' s are the coefficients 

of Fourier expansions in k^z for p' 2l and p^.The latter quantity is essentially 
the induced polarization at ^ , consequently the emitted power, AI , may be 
expressed as a sum of x(n)' s (n even), [Eq. (30b ) ] .For all physically 
possible (i. e. finite) values of <f> (n) AI must remain finite, hence the x(n)' s 
may not diverge. [It was shown in Ref. 21 that the <Mn)'s themselves 
strongly converge to 0 as n-»"». See also Table I. ] 


This requirement on x(n) places a restriction on the coefficient 
of v(n) in Eq. (Al8c) . Note first that in Eq. (Al8c) both summation 
terms (including multiplication by u (n) or v(n), respectively) con- 
verge rapidly to zero for n large. Thus, since u(n) itself converges 
** 1/ni [Eq. (A5a)], the entire first bracketted term of Eq. (A13c) 
converges. But, as can be seen in Eq. (A5b) , v(n) diverges *■' (n-1)!, 
placing a restriction on its coefficient. In fact, the required 
boundary condition can only be. guaranteed by setting 


The remaining constant, C(l), may now be determined. Displaying 
the subscripts as in Eq. (Al), we have: 
from Eq. (Alb), 



(A19) 


x + (!)-x_ (1)+A(0) x (0 ) = cf>(0) 


(A20a) 


from Eq. (Al7b) 


x + (l) = C + (l)u + (l)+D + (l)v + (l) 


(A20b) 


x (1) = C (l)u_ (1)+D_ (l)v_ (1) 


(A20c ) 


and from Eqs.(A18a) and (A17), 


x (0 ) = C + (l)+D + (l) 


(A20d) 


x (0) = C (1) +D (1) . 


(A20e) 


These equations may be routinely combined to obtain 


<f> (0 ) + ? (-l) j <f<|> + (j)u + (j) - <f>_ (j)u_ (j)} 



x (0) 


(A21) 


Substitution of definitions (A2) gives the result of Eq. (37) of the text. 


Rate Equation Approximation : 

The Rate Equation Approximation (REA) for a two level system 
interacting with a standing wave field is that approximation in which spatial 
population variations arising from the interaction are neglected. In extending 
this approximation to a 3 -level system we ignore contributions of the 
polarization arising from these terms (i. e. x n = 0, n > 1). The difference 
equations [Eq. (A1 a-cj) considerably simplify and we have 


x ( 0 ) +A (l)x (1) 
+ + 


REA 

<MD ' 


(A22a) 


x + (1) -x_ (1) +A (0 ) x (0 ) 


(A22b) 


Solving these equations directly, we find 


REA 

x (0) = 


REA 
4>(0) + 


REA 

<f>- (1 ) 

A (1) 


REA 


+ 


AjlTT 


A(0)+ a^tt) - nri 


(A2 3) 


REA 

where the <J> (n) are obtained from Eq. (A2)and E< 3* (35) , using the 
RBA vanuaes, A^ble I, Eqs t (T20) - (T22) . 


-45- 


Appendix B: Solution of the Second Order Linear Homogeneous 

Difference Equation (A3) : 

Consider Eq. (A3) , 

x+ (n+1) -x^ (n-1) + A + (n)x + (n)=0 , n>l, (Bl) 


where 


A + (n) = + — L(+n) , 


L + (n)=1 


r Ai 2 ~ (ki+ n k 2 ) v +iY 2 i' n odd 
A l~ (ki± n k 2 ) v+ iYic' n even 


(B2 ) 


(B3 ) 


In second order homogeneous difference equations, such as 
Eq.(Bl), one may distinguish two classes of solutions, character- 
ized by the behavior of x H (n) for large values of n. In the 

present case the nature of these types of solutions is readily 
ascertained, noting that for large enough n (for v^O), 

L+ (n) + +nk 2 v, (B4 ) 

so that 

2ik_v 

. A + (n) - — ^ — n. (B5) 

— p o 

To obtain the first class of solutions, assume that for large n 

xj (n+1) »x? (n-1) . (B6 ) 
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Then 

xj (n+1) 


(n) 


2 ^ 2 v 


so that 

x + (n) « (n-1) 



(B7 ) 


(B8 ) 


Note that x + (n) for large n, in accord with initial assumption 
(B6) . We denote solutions of this type by v + (n). Because of 
their close connection with Neumann functions we refer to solu- 
tions of the latter type as N-type solutions . 

The second class of solutions is obtained by assuming 

xj(n-l) >>x” (n+1). (B9 ) 


Then 


X J (n) 


x+(n-l) 



(BIO) 


so that 

x + (n) « 1 / g ° \ n . 

n! \ 2ik 2 v j 


(Bll) 


Note that x + (n) ->0 for large n, in accord with initial assump- 
tion (B9). These convergent solutions are denoted by u + (n). 

Because of their close connection with Bessel functions we refer 

( 21 ) 

to solutions of this type as B-type solutions. 


To express u + (n) 

define the quantity 
u + (n) 


W + (n) 


u + (n-1) 


in a form convenient for computation we 


(B12) 
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Then from Eq. (Bl) we obtain an expression for W + (n) in contin- 
ued fraction form, 

± iS 0 /2L[ ± n] (B13) 

w ± (n) = 

B 2 / (4L[±n]L[±(n+l)]) 

1 2 

fi 2 / (4L[± (n+1) ]L[± (n+2) ] ) 


l±ie Q W ± (n+3) / 2L [ + (n+2 ) ] 


Then 

n 

u (n)=II W + (k) , 

k=l ' 

where we have set u + (0)=l without loss of generality. 


(B14 ) 


Special Case y 2 =Yq and A 2 =0 


For the case y 2 =Yq and , u + (n) may 

in terms of J (z) , the Bessel function 
T 

complex order t: 


u + (n) = 


1 J t(n) ( - 6 q/V> 

J T(0) ( -W> 


be written analytically 
of real argument z and 


(B15) 


where 


t(n)= +L(±n)/k 2 V. (B16) 

As before, Eq. (Bl 5 ) is normalized such that u + (0) = 1. 

Equation (B15) follows from Eq. (Bl) which, in this special 

case, reduces to the well-known recurrence relation for cylin- 

( 21 ) 


drical functions. 
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Appendix C: Algorithm for Computing Continued Fractions 


Let 


C= 


1 + a. 


(Cl) 


Define 

P 


n 


n 


1 + a. 


1 + a. 


1 + a. 


(C2 ) 


1 + a. 


1 + a, 


1 + a 


n 


One can then show that P and Q obey the following difference 

n n 

equations : 


P = P , + a P _ 
n n-1 n n-2. 


(C3a) 


Q_ = Q__ n + aQ_ 

n n 1 n n 2 , 


(C3b) 


with 


P_1 = 1,P Q = 0 


(C3c) 


Q -i = 0,Q o = 1 - 


(C3d) 


From the recurrence relations for P an< ^ Q one can readily 

n n 1 


3 > 
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compute P n and for any value of n. The computation is halted 


when 


(P n+l/W - ( VV 


< 10 ” 6 . 


(C4 ) 


We then set 


C = 


n+1 . 


(C5 ) 


'n+1 


1 . 


267 (1966). 
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Figure Captions 

Figure 1 (a) . Energy level scheme considered in this paper. 
Note that is variable and may be greater or less than 

unity. 

Figure 1 (b ) . Typical gain at for standing-wave field 
at ^ detuned from resonance. The broad and narrow dips in 
the laser transition gain profile will be treated in a separate 
publication. The dips in the coupled transition are treated 
here. The insert depicts the dios as they occur in the absence 
of the background ( <R 0 2 > t * ie text) . In the example shown, 
the background (<R^^>) is unaltered by the standing wave field. 

In other cases both terms may be influenced by it. 

Figure 2 . Possible experimental arrangements for observ- 
ing laser-induced line narrowing. 

(a) Transmission - a tunable probe field is coaxial 

with standing wave laser field £ 2 (^ 2 ^* A filter blocks E 2 . The 
output intensity of the probe field is monitored as a function 
of . 

(b) Spontaneous emission - the fluorescence spectrum 
emitted along the axis of the laser cavity is studied by means 
of a high resolution scanning Fabrv-Perot inteferometer . The 
filter transmits only . 

The lineshape effects will be observed in both arrangements. 

Figure 3 (a) . Plot of the spatially averaged population dis- 
tribution vs. axial velocity for the upper level of a two-level 
system interacting with an intense standing wave field on reso- 
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nance . Note the fine structure appearing in the central dip over 
a narrow range of velocities about v=0. 

(b) . Typical curve of gain at as a function of the 

detuning from the atomic center frequency of the coupled transi- 
tion for the same casfe as Figure 3 {a) . Comparison indicates that 
structure similar to the population fine structure occurs at 
the wings of the narrow central dip of the profile. The central 
dip itself originates from complex interference effects dis- 
cussed in Section IV. 


Figure 4 . Gain at background omitted, for the case 

li) l~ t0 2 // ^ aS •'■ n ^^ uencec ^ ky an intense standing wave field 
detuned from resonance. Note that the dip on the right is split. 
For weakly saturating values of E^ r the splitting disappears. 

This new effect should not be confused with another effect which 
appears in the background term even for weakly saturating values 
of E 2 [see Figure 7(d)]. 

Figure 5 . Intense standing-wave field on resonance ^2=0) . 

(a) . < ^02 > versus for k^=k 9 . The additional structure 

on the sides of the central dip emerges as I increases. The 
latter structure is closely correlated with the structure in the 
space-averaged population distribution (Figure 3a) . 

(b) • <R oi > versus for the same conditions as in (a) . 

Note that this structure would be present even if the medium were 
transparent (n 2 =n Q ). 

(c) . <r q 2 > versus A i for k 1 =k 2 /2. The splitting observed 

in the central dip also occurs when E 2 is off resonance [Figure 
6(c)] and is a traveling wave effect. 


<- 


P-3 


(d) . <r qi > versus A^ for same conditions as in (c) . 

This structure also occurs in the detuned case [Figure 6 (d) ] and 
is a traveling wave effect. 

(e) . <R 02 > versus for k^=2k 2 . The complex structure 

observed within the broad central dip is the result of interfer- 
ence effects due to the standing wave nature of and does not 
appear when is detuned [Figure 6(e)]. 

(f ) . <r qi > versus A^ for the same conditions as in (e) . 

The structure here is also the result of standing wave interference 
effects and does not appear when E 2 is detuned [Figure 6 (f ) ] . 

Figure 6 . Intense standing-wave field well detuned from 
resonance (A 2 /y 02 =7.5k 2 /k 1 ) . 

(a) . <r q 2 > versus A 2 . for k^=k 2 » Two Lorentzian resonances 

occur at Q + =0 of FWHM r + , respectively, in accord with the analytic 
expressions Eq. (50c) . 

(b) . <R oi > versus for same conditions as in (a) . Note 
that E 2 has negligible effect on the background. 

(c) . <r q 2 > versus for k^=k 2 /2. A Lorentzian of width F + 

occurs at ft + =0. The resonance at ft_=0 splits into two, with increas- 
ing splitting for increasing values of I. 

(d) . <r qi > versus A^ for the same conditions as (c) . The 
background curve exhibits non-Lorentzian structure at £2_=0. No 
structure appears at £2 + =0. 

(e) . <r q 2 > versus A^ for k^=2k 2 » Lorentzian resonances of 

widths r + occur at S2 + =0 . 

(f ) . <r qi > versus A^ for the same conditions as in (e) . 

The flat curve indicates that E 2 has a negligible effect on 
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<R q . > when the standing wave field is detuned. 

Figure 7 . Weakly-saturating standing wave field. 

(a) . <r q 2 > versus A ]_ f° r ^2 =< ^‘ In t ^ ie weak li m it the 

resonances at A-^=0 are essentially the siim of two Lorentzians of 
widths r + and T_. 

(b) . <R oi > versus for A 2 = Note that for k^k.^, the 

background curve is unaffected by E 2 . However, for k^<k 2 , a 

non-Lorentzian structure appears about A 1 =0 . 

(c) . <r q 2 > versus for A 2 =7 ' 5k 2^ k l* The tW ° resonances at 

£2 + =0 are of widths T + . 

(d) . <r qi > versus ^ or A 2 = 7 * ^ k 2 //k l ’ T ^ ese curves are the 

same as the background curves of (b) except that the non-Lorentzian 
resonance at ft_=0 in the k 1 =k ? /2 curve is shifted. 

Figure 8 . Plots of <R Q2 > ( a ) and <R 01 > versus A^ for 

1=16, k^=k 2 , for the case in which E 2 is near resonance, A 2 /y q 2 =3.0 
These curves show features common to both the off-resonance and 
on-resonance conditions, indicating the complexity of the line- 
shape resulting from the interference effects treated in this paper 

Figure 9 . Comparison of <Rq 2 > versus A^ for 1=16, k^=k 2 for 
the exact solution [Eq. 43(b)], the IFA [Eq. (45)] and the REA 
[Eq. (A23)]. Note that structure in the exact curve on the wings 
of the central dip coincides with the residual structure in the 
IFA arising because of the fine structure in the population dis- 
tribution of level 2, although the structure in the exact curve is 
significantly enhanced. This enhancement results from spatial 
interference effects and does not appear in the REA curve. The 
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narrow central dip in the exact curve is in part due to the field 
interactions between and E 2 in the absence of spatial variations 
in the population of level 2, as can be seen by comparison with 
the central dip of the REA curve* The exact curve, however, is 
further influenced by spatial interference effects. 
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